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A su mmar y is given of the principal results thus far obtained 
from studies of the nonsteady motion of a lifting surface in an 
incompressible fluid; the methods followed by various investigators 
are Indicated. 

The aerodynamic problem of the nonsteady motion of an airfoil 
has been the subject of numerous investigations, which in 20 years 
have brought a degree of development such that an entire branch 
(that of two-dimensional motion) may be said to have been completely 
solved. The mass of existing publications is very large »nii a mon g 
these naturally many overlap; moreover, the study of the problem has 
produced a variety of methods such that the same phenomenon is 
endowed with rather diverse aspects that, although of considerable 
speculative interest, do not always facilitate the task of those who 
wish to learn only the results of the research. A synthetic deduc- 
tion of the results thus far obtained is given herein^- in order to 
expound the principles of the various methods of investigation and 
particularly to collect the latest results in a form that is best 
suited to application. 

As previously stated, the two-dimensional problem has been 
exhaustively studied. In this field, the problems of unsteady motion 
reduce to computations that at times are laborious, but which can 
always be conducted without uncertainties of the approxi nations and 
which are considerably facilitated by existing tabulations. Only a 
few attempts have been made to develop a rigorous computation for 
the wing of finite aspect ratio inasmuch as the existing methods 
present approximate solutions that contain many inaccuracies. 

In part I of this report, the results relative to the wing of 
infinite aspect ratio are described, always with the assumption of a 
perfect and incompressible fluid with regard to the components of 


. *"Lo Stato Attuale delle Ricerche sul Moto Instazionario di Una 
Superficie Portante." Estratto Da "L'Aerotecnica" Vol. XXI, N. 9-10, 
Settembre-Ottobre 1941, XIX. 

-*-The analytic developments have been shortened and the rigor of 
the demonstrations has at times been, relaxed. 
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of the normal force on. the wing plane . In part II* procedures are 
indicated for the computation of the wing of finite aspect ratio; 
a comparison is made of the various approximations thus far employed. 
In part III, the results are given relative to the force component 
in the direction of motion (propulsive or drag). In part IV, the 
results of experimental investigations are considered. In each part 
an evaluation is given of the existing publications, which may serve 
as a guide to a more detailed study of any aspect of the problem. 

Once a reference system is established with respect to which 
the fluid is motionless at infinity, it is assumed that a fixed 
plane exists from which the points of the lifting surface are at 
distances that may always be considered very hwm. 11 with respect to 
the dimensions of the projection of the surface on the plane. The 
thickness and the curvature of the profile sire therefore considered 
infinitesimal, as are the displacements along the normal to the 
fixed plane. Almost all the investigations up to the present time 
on the nonsteady motion of a wing use simplifications that are 
derived from the preceding assumption and thus assume that the sin- 
gularities in which the body in motion is schematized are permanently 
contained in the fundamental plane. In this plane lie the orthogonal 
axes x and y, which cure displaced with respect to the fixed refer- 
ence and remain parallel to themselves with a velocity V par all el 
to the x^ixis but opposite in direction; the z-axis is at right angles 
to the x,y-axes. With respect to these axes, which follow the motion 
the wing, the relative velocities of the points of the wing are 
small with respect to V. It is thus assumed that the perturbations 
produced by the motion of the wing are sufficiently snail. Hence, in 
the relations that are used, all terms of higher order than the lin ear 
terms in the velocities induced by the motion of the obstacle in the 
surrounding fluid are neglected. This fundamental simplification 
gives the problem under consideration the advantages of the Hna»> 
theories (the most important of which is the principle of superpo- 
sition) from the effects of which by the analysis of particularly 
simple motions, solutions can be obtained that through Himr combi- 
nations make possible the study of motions of a more complicated 
character. 

The velocity V is assumed to be small compared with the velo- 
city of sound although investigations have been conducted that consider 
the compressibility of the medium (reference l). 


The variations in time or space must frequently be measured. 

The symbols occurring in the derivations herein are therefore defined 
in the following paragraph: 
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Let q fee any magnitude that is a function of a point and that 
for every point varies in time* The symbol dq/dx (or dq/dy) 
denotes the dependence of the lav of variation of q as a function 
of x (or of y), measured fey giving to q, at each point, the 
values that correspond to a certain fixed time t. The symbol 
dq/dt indicates the rate of variation of q at a fixed point 
relative to the x, y, and z axes, vhich in the absolute refer- 
ence system is displaced vith the velocity V. In general, with 
the assumed linearization, the derivative also represents the rate 
of variation of q at a fixed point of the ving because the velocity 
of a point of the ving vith respect to the moving axes is very small; 
if the gradient of q is not too great, the variation of q in time 
vill fee the same vhether measured at the fixed point relative to the 
moving axes or measured at the point ^lat follows the lifting surface. 
In order to express the rate of change in time of a q always 
measured at the same point of the stationary reference, the symbol 
d'q/dt is used. Because of the assumed linearization, the previous 
derivative coincides with the derivative that measures the change in 
q with time for the same material point. Inasmuch as the absolute 
velocities of the fluid particles are very small, if the gradient of 
q is not too large, the change of q with time vill fee the same 
whether measured at the fixed point or following the molecule. It is 
preferred in this report to use a distinctive sign in d*q to differ- 
entiate the local derivative from the derivative of the quantities 
that depend for the problem under consideration only on the time param- 
eter; for this problem the notation dq/dt is vised. 


PART I. - TWO-DIMENSIONAL PROBLEM 

1. The results of modern research on the two-dimensional problem 
will first fee described. The simplification Introduced fey the assump- 
tion that the phenomenon develops in the plane of the x- end z-axes 
is such that it can fee stated that each problem within this range 
can fee reduced to the computation of integrals that, with graphical 
procedures aided fey analytical considerations, can fee computed with 
suitable accuracy (limited to the approximation of the existing 
tables), vhich eliminates the expansion into series. In the three- 
dimensional case (the ving of finite span), except for some parti - 
ular problems, only approximate solutions exist. 

In order to simplify the expressions of the two-dimensional 
problem, the semichord of the ving is assumed to be of unit length. 

The abscissas of the leading and trailing edge of the ving are assumed 
to fee given fey x * -1 and x ■ 1, respectively. In order to define 
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the points of the profile, the parameter £ is introduced. The 
variation with x is shown by the relation 


X * - cos $ 


The values 0 and n of the parameter therefore correspond to 
the leading and trailing edges of the wing, respectively. 

The coordinate z, normal to x, is considered positive in the 
downward direction (fig. 1). The vertical component w of the velo- 
city of the fluid is therefore considered positive if turned in the 
direction of positive z. The difference in. pressure p between the 
two faces of the wing has the positive sign directed upward. The 
same sign convention is true for the vertical force P, the lift of 
a segment of unit chord. The moment M on the same segment is con- 
sidered positive if it is a diving moment. 

2. Condition of tangency. - The condition that the profile be 
impenetrable to the fluid is expressed by malting the relative velo- 
city of the fluid, with respect to the wing, tangent to the profile; 
or, in other vords, the absolute velocities V c of a point on the 
contour and Vf of the fluid particle in contact with the contour 
have the same projection on the normal n (fig* 2). Tnft«m»ch as 
only the linear terms in the coordinate z of the points of the con- 
tour and their derivatives, or in the velocities (perturbations) 
created in the fluid by the motion of the wing are considered, the 
component an n of the velocity may be supposed equal to the 

component w of the velocity of the fluid particle, which is parallel 
to z. By using the simplification of neglecting the quadratic terms 
in the computation of the projection on n of the velocity V c , 

which has the components V and dz/dt, there is obtained 


w 



dz 

3t 


( 1 ) 


3. Circulation and pressure on profile in steady motion. - Under 
the conditions of steady motion, it is known that the relative velo- 
cities of the fluid with respect to the profile, even if the profile 
is considered to be of infinitesimal thickness, are, in general, 
different on the two surfaces. Xn the motion under consideration, 
if Vi and V2 are the velocities at the corresponding points of 
the two surfaces (fig. 1), the pressure rise p between the surfaces 
is given by the Bernoulli equation 
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1 y_+y- 

P - 2P(T 2 2 - t2 i> - P(V2-Vl) ~~ 2 ~ - PVu 

where it may "be assumed, because of the linearization hypothesis, that 
the average of the velocities Vi and is Y, and u denotes 

the difference between them. 

If the wing is in the positive aspect, the smallest velocity is 
found at the bottom surface. The velocity difference u is con- 
sidered positive under these conditions and the positive pressure p 
is therefore in the opposite direction to positive z. In order to 
represent the field created by the wing, the skeleton of the wing 
composed of a vortex film is considered; this system of singularities 
is capable of giving the existing velocity increment u between the 
two surfaces if the circulation in an element dx is equal to u dx. 

4. Circulation and pressure on profile in unsteady, motion. - 
For the case where the motion is unsteady, the velocity increment may 
be represented by a vortex distribution of intensity g - u. The only 
difference, when compared with the preceding case, is in the fact tha t 
the discontinuity in the velocity field exists not only on the points 
of the wing but also in the wake behind the wing. In order to describe 
this phenomenon, consider two fluid layers that pass above and below the 
profile. Because of the linearization assumption, the distances of the 
points of the profile from the x-axis can he neglected in the following 
discussion. It can easily be verified that the terms of the second 
order will therefore be neglected. The difference between the momen- 
tums of the two layers is computed (fig. 3) 2 at the time t 0 and at 

the time t 0 + dt and therefore the variation that the difference has 

undergone in the interval dt is also computed. This variation, 
divided by the interval dt, must be equal to the difference between 
the forces along x that are applied to the two layers. T "a «imi'-h as 
tangential actions do not exist (perfect fluid), and because no pressure 
difference exists on the anterior face (upstream of the profile), the 
previously mentioned momentums will be given by (pi-p2) dy ■ p dy. 

The difference in the momentums of the molecules of the two layers, 
which is given by the product of the mass and the difference in the 
velocities u, varies in the interval dt because: 


2 

The two layers are indicated in figure 3 by hatched lines in the 
two directions; the position of the leading edge of the profile is 
indicated by a semicircle. 
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1. New molecules come in contact with the profile^ and there- 
fore acquire the velocity increment u = g. 

2. The increment for a given point of the profile varies with 

time . 


If the graphs of g relative to t = to and t ■ tQ + dt on 

the same position of the profile (fig. 3) are plotted corresponding 
to the two causes, respectively, it is found that the area (which 
gives the difference in the momentums when multiplied by P dy) varies 
in the interval dt by the amount 




gV dt + dt 


if 


-1 


The first term is indicated in figure 3 by the obliquely hatched 
area (which, except for infinitesimals of higher order, is equivalent 
to a rectangle of base V dt and altitude equal to the value of g 
corresponding to the abscissa x of point P) ; the second tern is 
indicated by the vertically hatched area. Equating the impulse to 
the variation of the momentum and dividing by p dy dt yield 4 



( 2 ) 


The first term on the right side measures the pressure at P, 
which is obtained under the conditions of steady motion; this pres- 
sure depends on the local value of the velocity increment. By the 
effect of the other term, the pressure under the conditions of non- 
steady motion depends on the variation that g undergoes at the 
instant considered in the entire strip ahead of P. 


"The difference of the momentums of the two layers is observed to 
be zero before arriving at the wing because the velocity increment in 
front of the wing is zero if there are no other lifting surfaces. 

^This relation is obtained by another method in differential 
form in reference 2. 
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5. Vorticity of -wake. - If at the instant t = to the point P 
coincides with the trailing edge of the wing, equation (2) still 
holds, hut the first term becomes zero because a difference in pres- 
sure cannot exist where there is no wing surface. It follows that 

the increase in velocity g_ existing in the points of the wake, 

s 

which is zero in the case in which the distribution of g on the 
wing is independent of the time, is given by 

g s = " v dt 


where 

P 1 

K = J g dx 


is the circulation about the wing and the derivative is measured at 
the instant at which the trailing edge passes through the point of 
the wake under consideration. In other words, in the distance V dt 
that the wing moved in the time interval dt a vorticity is distri- 
buted equal to the variation that the circulation about the profile 
has undergone in the same time. This conclusion can also be derived 
from the principle of the conservation of vorticity. 

In figure 4, the diagram of the vorticities on the profile and 
in the wake is given for the case in which the wing executes a trans- 
latory oscillation with frequency Q related to L and V by 


C3L = 27 


The curve a refers to the instant in which the wing crosses 
the middle position; b refers to the position at the end. The 
scale of the vortices is indicated by assuming the vertical ma^mum 
velocity of the profile to be unity. For comparison, figure 4 shows 
the graph of the circulation g' corresponding to the maximum veloc- 
ity in the case of steady motion. The increment of velocity in the 
wake is considerable. In the theories of a wing in nonsteady motion, 
it is generally assumed that this velocity increment remains in the 
position in which it originated, or, in other words, that the vortices 
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shed by the wing maintain their original position and intensity 
unaltered in time. It is evident that the error inherent in this 
assumption, from which the real phenomenon certainly deviates, 
impairs to some extent the results of the theory. 

6. Bound and free vortices. - An original physical interpre- 
tation of equation (2) was given by Bimbaum (reference 3). By 
setting 


7 


P_ 

py 


( 4 ) 


€ 


Px 


1 

V 



(5) 


Equation (2) is written as 


g = 7+ e 


( 6 ) 


With the aid of this relation, the value of the total vorticity 
g at a point on the profile at which the velocity increment between 
the corresponding points of the top and bottom surfaces is divided 
into two parts : 

(1) The bound vorticity y, which is sustained by the aerody- 

namic action 

(2) The free vorticity e, which trails in the fluid in its 

relative motion and which therefore gives no pressure 
rise 

In order to clarify the relation between the free vortices and the 
bound vortices, equation (2) is differentiated with respect to x 
after g is expressed in terms of equation (6). Thus, 





be _ 

<5x 


d'C 

dt 


(?) 
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This relation indicates that the variation undergone in an inter- 
val of time h.y the vortex C at a point of the fluid is equal (and of 
opposite sign) to the variation undergone in the same time hy the 
bound vortex at the point of the profile in contact with it. In 
other words, the bound vortices leave at each instant, in the fluid 
with which they come in contact, an effect represented by a vorticity 
of intensity equal to that which they have produced at that instant. 
According to this representation, each of the inducing elements at 
the wing is considered in isolation and therefore produces a pres- 
sure rise expressed by equation (4) and has a proper vortex wake. 

The free circulation at any point on the wing or downstream of it 
is given by the sum of the circulations of the wakes corresponding 

to the bound vortices that are upstream of the point considered. 

The decomposition of the total circulation in the two vortex 
systems previously described is shown in figure 4. It can be seen 
that the graph of the free circulations on the chord are joined in a 

continuous manner with the system of the wake at the point where the 

free circulation is equal to the total circulation. 

7. Relations between circulations and normal velocities. - 
Ordinarily, the law of variation of w along the chord and with time 
is known from equation (1), with the aid of which these quantities 
are derived from the characteristics of the motion. The vorticities 
on the wing and in the wake must induce at each instant on the points 
of the chord the assigned w; that is. 


v(x) * s/ 

The integration is extended from the leading edge of the wing 
to the entire part of the wake in which the passage of the wing has 
created the vorticity discontinuity. 

It is known that, for the wing in steady motion, the field of 
motion would not be determined if the point of separation (trailing 
edge of the wing) were not fixed. This condition is also assumed 
for the wing in unsteady motion and is translated into the analytical 
condition that g be finite at the trailing edge of the wing. 

# 

The vorticity in the wake is connected with that at the wing 
by equation (3)'. When equation (8), which is completed by the condi- 
tion of separation and by equation (3), is solved the total vorticity 


g(x* ) dx 1 
x' -x 


( 8 ) 
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g is obtained. The pressures are then obtained by making use of 
equation (2). As will be shown in more detail in section 17, this 
scheme is followed by various procedures of solution by Theodorsen, 
who makes use of conformal transformations (reference 4), Schwarz 
(reference 5), who makes use of Betz's solution of equation (8), and 
Sohngen (reference 6). Wagner and Glauert also refer to the total 
vorticity, but determine the total action on the profile and not the 
pressure distribution (as do von Kalman and Sears, reference 7). 

When a different method is used, the free circulation on the 
chord and in the wake can be expressed as a function of y by 
making use of the integrated equation (7 ) and thus transforming 
equation (8) so that only the unknown function y appears in it. 
Then the pressures can be directly obtained by means of equation (4). 
This procedure was followed by Birnbaum (reference 3), by Kussner 
(reference 8), and by Cicala (reference 9). 

8. Acceleration potential. - A different interpretation of the 
same problem can be made on the basis of the acceleration potential. 
For a perfect and incompressible fluid, the equation of Euler 


a = - grad £ 

expresses the equality between the force of inertia, which corresponds 
to the acceleration a of the fluid particle, and the resultant of 
the forces that are transmitted to the particle by the medium surround- 
ing the particle. The Euler equation also permits stating that the 
components of the acceleration can be obtained from the function p/p 
through the same operations of differentiation with which the com- 
ponents of the velocity are deduced from the corresponding potential. 
The generating singularities are arranged on the surface of discontinu- 
ity (wing + wake) for the acceleration field whose potential satisfies 
analytical properties similar to those of the velocity potential in 
the same manner as for the velocity field. The singularities are 
arranged where the discontinuity exists in the pressures; that is, 
on the lifting surface. The simplification that is introduced by 
this concept is not, however, as great as might appear from the fact 
that the singularities (which are called the dipoles of the pres- 
sure) are limited only to the wing. In each case, it is necessary 
to pass from the acceleration field to the velocity field and this 
passage requires an integration through which the effects of all the 
preceding states of motion are felt and which, according to the 
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vortical representation, leaves a trace in the system of free vortices . 
In other words, the velocity field of the pressure dipoles depends 
on the history of the formation of the dipole. 

The relation "between the vortex representation of the phenomenon 
and that based on the pressure dipoles is already implied in the con- 
cept of Birnbaum of the vortices shed from all the points of the 
airfoil that is the only source of the bound vorticity 7 , which is 
proportional to the difference in pressure between the two faces of 
the profile. A method will be shown, on the basis of the concept of 
vortices, for the derivation of some fundamental relations that in 
other publications are justified with the procedure of the accelera- 
tion potential. 

The system of bound vortices induces at a point 0 at the 
instant t = 7 a velocity that is denoted T^; at a succeeding 

instant t = 7 + A T , at a point that occupies, with respect to the 
wing, the same position that 0 occu]$Led in the first condition, 
there will be a velocity differing from Y^ by an amount denoted 

by AV fl . The velocity AV n will depend on the variation in the con- 
stitution of the induction system or, with changed sign, will be the 
velocity induced by the system of vortices that, at the instant con- 
sidered, are freed from the bound system and constitute the trace 
that the bound system has left in the fluid with which it has been 
in contact. The quantity AV /at , as A7-K) as a limit, represents 

the derivative d v^ /dt. Hence, if it is desired to express analyti- 
cally the property that the actual velocity Yf °f the fluid parti- 
cles results from the sum of the velocities induced by the actual 
configuration of the bound vortices and from those freed at the pre- 
ceding instants, then 


ft 

dv , x 

* ‘ - 3t dT (9) 

J - Od 

where T a and Vf are computed at tne point and at the instant t, 
dv/dt being computed at the point of the fixed space being considered, 
and the system of. bound vortices is in the corresponding position at 
the instant 7 preceding t.. In the case of translational motion 
or, more precisely, in the linearization assumed, the terms of the 
second order in the velocities of the points of the wing relative to 
a system that is displaced with translational motion with velocity V 
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are neglected; the term dv a /dt is measured by having the point 
of induction at the abscissa x-V(t-T), or, with respect to the 
point 0 of the abscissa x to which the velocity Vf refers, 

in a position upstream, advanced by the quantity V(t-T) if V is 
assumed constant. 

The following relation can be obtained from equation (9): 


d'Vf „ ^la v 

f s- T 


( 10 ) 


The quantity V a is the velocity induced at 0 at the constant 
t by the bound -vortex film layer. From the instant t-dt to the 
instant t, the vortex system changes only in the number of the bound 
vortices leaving a trace and that occupy the position borresponding 
to the time t. The variation tha| the velocity Vf of the fluid 
has undergone at the point 0 of the fixed space because of the 
effect of this change (local variation d'V^/dt dt) is that which 

would be measured by leaving the vortices y in position and moving 
the point of induction of the segment dx = Vdt so that the point 
of induction occupies the position that corresponds to the instant t; 
the variation that is measured in this manner is expressed by 
dVa/dx dx and is thus added to equation (10), which equates the two 
expressions of the variation. 

Relations (9) and (10) are equivalent 5 to the relations derived 
by Possio, which are based on the concept of the acceleration poten- 
tial, if it is •assumed that the stationary field C, defined in 
reference 10, coincides with the field produced by the system of 
bound vortices in a uniform stream V. Inasmuch as the concept of 

^Equation (10) is derived directly if the acceleration field is 
considered to depend only on the actual values of the pressures on the 
profile and, with the pressures equal, has the same configuration as 
though the motion were steady. In this case, according to equa- 
tion (4), the actual pressures would be obtained by having only the 
actual values of y on the profile: the velocity of the fluid would 

be V a and the acceleration that, by the assumed linearization, is 
computed from d'V a /dt (and not from the derivative formed in following 
the fluid particle) would be expressed by d'V a /dt = dv a /dx Y (without 
dVa/ot because this virtual field is steady). When the acceleration 
of this virtual field is equated to the effective acceleration, which 
is written d'Yf/dt, equation (10) is obtained, which when integrated, 
gives equation (9). 
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the acceleration potential lends itself less to physical interpre- 
tations than the concept founded on vortices although some simpli- 
fications can "be obtained in the analytical development* it is 
preferred, in the following sections, to employ the classical method 
of description. 

9. Cases of total zero circulation. - If the total circulation 
about the wing is constant in time and the system has no vortices in 
the wake, the induced velocities can be computed on the basis of the 
total circulation about the wing at the instant considered in the 
same mann er as for steady motion. From the analysis relative to this 
case (reference 11, p. 185), it is known that if the circulations are 
represented by one of the following functions 


$ 

g-L = cot g - 2 sin d 
g n = - 2 sin n£ (n 



(ID 


the corresponding velocities are respectively given by 


wj_ = 2 + cos $ 


w n = cos n£ 


( 12 ) 


Inasmuch as the values of g n represented by equation (11) 
satisfy the condition 



the preceding result holds for any motion because no wake exists down- 
stream of the body. Thus, 


w -lA n w n 
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where A n are quantities independent of d hut functions of time. 

The total circulations are represented hy the corresponding sum 

g = 2 A n g n (13) 

The corresponding hound vortices are represented hy S y n} where, 
as derived from equations (5) and (6), 



In general, the values of w can he expressed hy expanding in 
a Fourier series: 


w = | A 0 + 2A n cos n£ (16) 

The coefficients A n , general functions of time, can he obtained hy 
harmonic analysis hy setting 


w = w - ZA n w n 


There is easily obtained 


w 


= | (Aq - Ai) 


(17) 


It is therefore concluded that, if the vertical velocities on 
the profile are developed in Fourier series, the circulations and 
the pressures can he directly computed with the aid of equations (13) 
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and (14) if Aq = A]_. In the general case, it will he necessary to 
sum the circulations (equation (13)), which correspond to the 
velocity 

w = — I w(l - cos d) dd (17* ) 

n Jo 

which is a function of time, hut is constant over the chord. 

10. Case of velocity w constant over chord. - Applying equa- 
tion (10) for the potential on the z-axis, inasmuch as w is the 
projection of V^, yields 


d 1 w _ dVa 
V dt 3x” 


(18) 


where Va is the projection of the vector on z. 

Because at each instant w assumes, for all the points of the 
chord, the value w, the following relation is obtained: 


d 1 w _ dw 
dt dt 

and, therefore, from equation (18), integrating along the chord yields 


dw 
V dt 


x + 


C 


where C is a constant with respect to x, hut is a function of time. 
From the analysis of steady motion, it is known that the distribution 
of the vorticity capable of giving velocities satisfying the condition 
at the edge is given by 


7q = 2 sin d + 2 C cot - 


(19) 
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In order to Investigate the dependence of C on the condition 
of motion, the elementary case is considered in which w, which has 
been zero for an indeterminate time, suddenly acquires, at the 
instant t' , the value Aw and maintains it unchanged. For such 
a case, for t > t' , C, which must be proportional to Aw, may be 
put in the form 


C = (1-R) Aw (20) 

S 

where R is a function of the space passed through by the wing 
since the instant t* . Because, with the passing of time, the phe- 
nomenon tends toward the steady conditions for which the circula- 
tion tends to assume the distribution 


_ ■d 
7 0 = 2 Aw cot g 


( 21 ) 


the asymptotic value of R must be zero. The law of the variation 
of R was studied in one of the first publications on the wing in 
unsteady motion (reference 12, which gives a resume of the work of 
Pistolesi to whom reference is made). The more general case is 
obtained from the elementary case by superposition of the effects. 
According to equation (20), the second term of expression (19) is 
decomposed into the asymptotic term (equation (21)) and a term that 
contains the function R and represents the distribution of the 
circulation that would be realized if w possessed, from an indeter- 
minate time, the value Aw up to the instant t 1 , and then for 
t > t' , w = 0. The pressures corresponding to this term, which 
represents the effect of the preceding variation and diminishes to 
zero with time, can be referred to as "transitory pressure." 

In the computation of C in a general case, summing the effects 
of all the increments that w has received from the start of the 
motion, when it may be supposed w = 0 (so that 2 Aw = w), there is 
obtained 



R dw 


( 22 ) 
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The integral is taken as the sum of the products of all the vari- 
ations that w has undergone in the preceding instants (the procedure 
is not modified whether the variations are abrupt or gradual) for the 
corresponding values of R. The distance traveled by the wing, 
measured from an arbitrary origin, is denoted by s; the value of s 
for the position at which the pressures are measured is denoted by 
s 0 and the value of s for the distance s 0 -s referred to the semi- 

chord is denoted by o • The values of the function R(o) are given in 
table I. These values have been obtained from the recent tabulations 
of Kussner and Schwarz (reference 13). 

11. General solution. - When the expressions of the two pre- 
ceding sections are collected, the resultant pressures for a general 
case of motion may be computed. If the values of w are expressed 
by the series of equation (16), which is put in the form 


v - v + EA n w n 


(16a) 


the pressures corresponding to the first term, according to equa- 
tions (19) and (22), are given by 



whereas, the values of w that constitute the summation of equa- 
tion (16a) correspond to the pressures 


p * PV Zr n 

where y n is expressed by equation (14). 

It is convenient to divide the total pressures into a part that 
depends on the history of the motion and has been denoted as the 
transitory pressure given by 


_P_ 

PV 


£ 

2 cot ^ 


R dw 


(23) 
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p, -which depends only on the actual -values of the 
characterize the motion and is denoted hy the 
pressure, given by 

p _ £ dw 

py = 2 V cot 2 + 2 sin y + Z 7n 

When equations (14) and (17) are used, this relation assumes 
the form 


and into a part 
parameters that 
•* instantaneous" 


£ 

PV 




A’ 


n+1 


2n 


^'n-l ^ 


sin n 




(24) 


For the computation of the pressures on the profile, the fol- 
lowing operations are therefore required: 

(I) From the law of motion there is obtained, with the aid 
of equation (1), the expression for w as a function 
of time and of the coordinate 

(II) By developing w in a Fourier series, the values of A 
as functions of time are computed. 

(Ill) From the instantaneous values of A and the derivatives, 
the values of ?S are computed from equation (24). 

When v has been determined from equation (17), the 
values of p are computed by equation (23) and hence 
the resultant pressure p » p + f . The total force 
and moment are obtained by the simple integrals 



^1 

M = px dx 

J-l 



NACA TM 1277 


19 


The formulas given here are also valid without change if V 
varies from one instant to the next. In such a case, it is con- 
venient to assume as a reference variable instead of time the dis- 
tance s traveled by the wing. Therefore 


A' n = ; where ds = V dt. 

In the calculation of the transitory pressures, the expres- 
sion (23) is usually computed by graphic integration. When the dia- 
gram of w as a fraction of 0 is known and the value of R is 
obtained from table I, it is convenient to draw the graph of w(R) 

(fig. 5), which is obtained by laying off, for each ordinate w, 
the corresponding value of R for the same position. The area 
enclosed by this curve, by the vertical axis, and by the two hori- 
zontal lines through the ends (one of which is the axis w = 0) 
represents the integral of equation (23) (cross-hatched area in 
fig. 5). If when tracing around the contour from a to b the 
area is on the right, the area is considered positive; in figure 5 
the value would be considered negative. Even though the graph of w 
presents abrupt variations (as in the case of the figure), no compli- 
cations are thereby introduced. 

When it is possible to proceed by the analytic method in com- 
puting the preceding integral, it is convenient to assume as the 
variable of integration the distance s. Equation (23) is there- 
fore written as 

s Q _ . 

||r(s 0 -s) ds (23*) 


or as 


i - 2 cot | f“-J 5 B(O) d0 (23") 

12. Example of application. - Let it be assumed that the wing 
always displaced with constant velocity in magnitude and direction 
undergoes a sudden small rotation (of small amplitude) about a point 


JL 

PV 


- 2 cot g 
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of the chord. Let there he determined the law of the variation of 
lift of the focal moment at the instant of the rotation. 

It is convenient first to assume that the rotation occurs in 
a finite interval of time and then passes to the limit to let the 
interval, approach zero. Let a he the angle of rotation^ cp the 
final value of a, and $ = do the coordinate of the axis of rotation. 

While the wing travels through the distance from s « 0 to s = A, 
the angle a increases continuously. Let 


for 

s 

<0, 

V 

IB 0 

for 

0 

< s < A, 

V 

= aV + (cos * 0 “ cos 

and for 

8 

>A, 

w 

- qpv 


Hence , 


AO 

2 


V^a 


da 

ds 


C08 



_ v ^ 

ds 


w 

* 


( 


Y a + r 


da\ 
ds / 


r 


| + cos 


for s between zero and A the rotation is w » V<P* In the first 
phase , 
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and in the second phase, dw/ds * 0. If s » sq (the position for 
which the pressures are measured), let Eq, E'q ... "be the values 
of E and the derivatives for o » s 0 . Equation (23* ) is then written 
"by expanding E in a power series in s: 


-P » 


2PV 2 




Then 


PA 


da . 
di ds 


0 


9 


(25) 


dfa 
ds 2 


ds = 0 


PA 


JO 


d 2 a . 

S jr dS 

ds^ 


9 


and moreover, the quantities 


PA 

0 


s“£ds 

ds 


'A o 

^ ds 
ds^ 


0 


for A approaching zero become zero at A n . 
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Thus, from equation (25), for A approaching zero, 

-P - 2pcpV 2 cot | (Rq - r R'o) 

From this relation it is concluded that the transitory pressure 
decreases according to the function R if r * 0 or if the rota- 
tion occurs about the neutral rear point. If r is different from 
zero, a term is added in the law of •variation of p that decreases 
as the derivative of R; this term corresponds to the pressure dis- 
tribution that is created on the wing in uniform rectilinear motion 
and that executes an instantaneous displacement in the direction 
normal to the trajectory and then continues with the initial speed 
and direction. 

The pressures ? after the rotation are given by 

f = cot 


Hence, the lift after the rotation is expressed by 


P 



(p + p) dx = jtp'CpLV 2 (1 - R + rR« ) 


(25d) 


where L is twice the chord of the profile and R and R» are the 
function of table I and its derivative, respectively, both approaching 
zero with an increase of the independent variable, which is repre- 
sented by the distance of the actual position from that at which the 
rotation has occurred. 

In take-off of the wing, that is, when the wing starts its 
motion from rest, 


w = Vcp 

In the interval 0 < s < A, V passes from zero to its final 
value V Q . For A — >0, 
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''A 

lim E ds g = R 0 <pV Q 


and the lift is therefore still expressed by equation (25a), in which 
the factor in parentheses reduces to 1-B (Wagner's case). For the 
moments, because the pressures are always proportional to cot d/2, 



In the preceding computation, the impulsive pressures that are 
generated at the instant of the rotation and cease when rotation 
has occurred were considered. The values are therefore immediately 
obtained with the aid of equation (24). 


13. Computation in finite terms of instantaneous pressures. - 
For the determination of the instantaneous pressures, a convenient 
expression is given by Sohngen, by means of which these expressions 
are obtained directly from w without expanding in a Fourier series. 
The expression, modified to conform with the notation used herein, 

may be given in the following form: 

o 


p/p V « A 0 cot | | sin d 


where 


H(^,s) ■ w(x*,s) x 3£ dx ( 27 > 


and where x = - cos d is the coordinate of the point to which p is 
referred and x' ■ - cos d* is the variable of integration. 


n* 


H( d* t B ) dd' 
COS d’ - COS d 


(26) 
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The lower limit of the integral in equation (27) is arbitrary, 
so that H is defined except for an additive constant. This arbi- 
trariness does not affect the results because 


n« 


\ j 


0 



- 0 


The proof of equation (26) can be given in a manner that is not, 
however, entirely satisfactory from the mathematical point of view 
by substituting the expression of equation (16) for w in equa- 
tions (26) and (27) and verifying that the relation thus obtained 
agrees with equation (24). 


The practical computation of the integral that occurs in equa- 
tion (26) presents difficulties for the sing ular ity of the function 
integrated at the point d = i>Q. For cases that are encountered in 

practice, however, by dividing the chord into a certain number of 
strips it is found that in each strip H can be represented by a 
combination of a few terms of the type cos nd. Hence, to simplify 
the applications, Sohngen gives the following formula, which in this 
case permits conducting the computation in closed form : 


P* 


sin d 


U 


0 


f dd' 

COS d* - COS d 


« 


_ cos nd ln 
2 


£ 1-COS (d-d2 )H ["l-COs(d+dl)"] 

[^(^^[r- T osid-d!)] + 


2 ^ ^ — sin(n-u) d (sin udg - sin hd^) + (d 2 - d^) B in nd 
u=l 


(28) 


where 


0 

for 

0 

<d< d, 

cos nd 

for 

♦l 

<d< d 2 

0 

for 

^2 

< d < « 



MCA TM 1277 


25 


In each case, hy means of a few terms of this type, the func- 
tion H can he expressed and the pressures can therefore he com- 
puted as a sum. The computation of the resulting actions involves 
easy integrations. In order to simplify this part of the computa- 
tion, Sohngen gives the expression 



-cos 


(&-a)2 sin nd d£ 


1 

—ip cos n 
n c 


£ 

n 


sin na + 


n-1 . 

| £ + i (oos no . oos n4) 

1 


The use of these expressions will he clear from the examples 
that follow. 

14. Examples of application. - Let the expressions of the pre- 
ceding section he applied to the determination of the Instantaneous 
pressures corresponding to the rotation of the elevator; that is, 
it is assumed that the forward part of the wing, corresponding to 
■values of $ between 0 and remains Immovable while the rear 

part rotates rigidly about the hinge located at the point of the 
abscissa x = - cos £q. The angle of the elevator is denoted by p 
(positive downward) and the primes denote the derivatives P' ■ dp/ds, 
P" = d 2 p/ds 2 . Hence, 


z * P(cos ^0 ” cos £ ) 


and therefore, from equation (l). 


w/v » p + p* (cos - cos $ ) 


H = w+V 






(P'+P" cos^q-P” cos sin£ dd «Bq+B^ cos £ + B 2 cos 2# 


vJ 
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where 


( IP” 

B 0 « V (P + 2P* cos $ o + 2 P" + T" 008 2»0 

B x - - V (2p« + P n cos $ 0 ) 

B 2 * “p- cos 2$ 0 


For each of the three terms of the expression of H, equation (28) 
is applied with n * 0, n * 1, and n * 2, respectively, and with 
^1 = ^Of ^2 = «• Making the substitutions yields 


pit 

. n 5 d* 

sin v CO g _ COB $ 


1 (B 0 + Bi cos $ + Bo cos 2$) In k Z 008 + *0) + 

2 1 - cos ( $ - * 0 ) 


(it - $ 0 ) (®l sin ^ + B 2 sin 2$) - 2 B2 sin £q sin $ 


On the other hand. 


-•■I * 


2( It - *q) O0t 

d$ = - (P ,+P* cos $ 0 ) + “p sin 
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There la thus finally obtained from equation (26) 


cot Z + 


« p/2 PV 2 =^0(rt-$o) * s '[ 8i " *° + <-»0> *o]j 

j^2p’ (it-^ 0 ) + P” (*-^q) 003 + 5 8in «b] 8in ^ " T* sin 2 ^ 

i r 1 jfl ‘ 1 -COb(i>+^o) 

||_3 + 20’ (cos $ 0 - cos $) + | 0” (cos * 0 - cos $ ) J In 1 _ cos ^_^ o y 


As a second example, the pressures in the case of the stationary 
gust are computed; that is, the wing is assumed displaced with the 
velocity, which is constant in magnitude and direction, encountering 
air layers that move in a vertical direction perpendicular to V 
with velocities that, at each point of the fixed space,. are main- 
tained constant in time. The graph of w along the wing trajectory 
(shape of the gust) is assumed given and P is a general point of 
the chord that is indicated by the positions corresponding to t = to 
and t = to + dt in figure 6. For the point considered, the value 
Ox 

of J w dx is represented at the first instant by the obliquely 

hatched area and in the succeeding instant by the same area increased 
by the horizontally hatched strip and decreased by the vertically 
hatched part. The intervening variation in the interval considered 
will be represented, except for infinitesimals of the higher order, 
by the quantity (w^-w) ds. Hence,- 


H « w + 



V 1 


The quantity H is therefore constant for all points of the 
chord. The same conclusions evidently hold for the point P if 
the point P has not yet entered the gust. (The only difference 
with respect to the preceding case is that in this case w = 0. ) 

Hence, from equation (26) (the expression for p), the only nonzero 
term. will be cot $/z. It is therefore concluded that for the wing 
that crosses a stationary gust, whatever the form of the gust, the 
pressures are distributed proportionally to cot $ /2. The same result 
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could also be deduced by considering equation (18) > for which, in 
this case, the first member should be zero from the hypothesis that 
w does not vary in time. The result given by Kussner that the 
aerodynamic actions on the wing that enters a stationary gust have 
a resultant passing through the focal point is thereby obtained. 

This result holds for the case where the values of w do not vary 
locally; in general, in agitated air the velocities vary rapidly 
with time. 

By following th3 analysis of the problem of the stationary gust, 
the case of the elementary gust (the step diagram in fig. 7 ) is first 
considered. During the time in which the front of the gust lies 
within the wing, . 




cos £ d£ 


2 3ln ^' Aw 
n 


v 


- sin 1 

it 


Aw 


where x*.= -cos is the abscissa of the gust front. When the 
entire wing is enveloped by the velocity Aw, 


A o = 2 Aw A]_ = 0 


The instantaneous pressures are therefore represented by 


— Q 

«p = 2pV Aw £ ' cot - in the first case, 

2 


and by 


£ 

|5 = 2pV Aw cot — in the second. 
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For the computation of equation (23), only during the phase of 
the crossing of the gust front is dw different from zero, and there 
it has the value 


dw = — (l-cos ) d^' = ^.JEEL dx' 
it n ** 1-x* 

If s denotes the distance of the midpoint of the wing at the 
actual position from the gust front, for which the values of p are 
measured, at this instant 


-it p/pv Aw cot % ~ ^ 

U 


1 R , 

R(s ” x> } dx> 


The pressure may therefore he computed as a function of R, in 
the case of the gust, by means of a simple integration. The resultant 
pressures may therefore be expressed in the form 


P 


• # 
p + p = PY Aw Rq cot g 


(29) 


where Ej is a function of the distance s^ ■ s '+ 1 of the front 

of the gust from the leading edge of the wing. The function is 
evaluated in table 2; the values are obtained from reference 13. 

When sj_ is negative, then evidently Ej_ = 0. 

From the solution relative to the elementary case, the solu- 
tion for a gust of any shape can be obtained by substituting in equa- 

J 'w 

R^ dw (taken as 

0 

the sum of the products of all the variations that the values of w 
undergoes for the corresponding values of R^) . 

15. Profile in harmonic oscillatory motion. - The same relations 
permit solving the case of harmonic motion. Assume V constant. 
Using the complex variable notation yields 
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where H is the frequency and V Is a function of d hut hot 
of t. Equation (26) for this case yields 


n p/apYe 1 ^ * cot 


» p* r* 

| W &■& - sin d 

Jo Jo 


W&jJ dd< 


COB V* - COS d 


r* 

oo sin d I 

Jo 


d*' 


cos d ' - cos 


-jr 


W(x) dx 


(30) 


where 


w * iy x* = - cos ■d' 


For v/ 


w * We*^ = We^ 8 (w constant) 


(For the position s - o , 


w(s-c) » wtsje" 41 ^ * We 1 ^ 1 


-ooo 

e 


(Hence, from equation (23?*), 


p/2PWe^ b cot | « a) R(o)e w O do (31) 

Vj 0 


The same problem can he attacked hy making use of the general 
relations initially given* By procedures that are developed in 
numerous publications and that, in part, are herein presented, rela- 
tions are arrived at that are equivalent to equations (30) and (31). 
In this manner, which is more rapid than direct computation, it is 
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found that the quantity that enters the second member of equa- 
tion (31), a function of the reduced frequency uT = Ql/ 2V, is iden- 
tified with the quantity denoted by \ in references 2 and 9 and is 
therein expressed by means of the 0ankel function of parameter u5 : 


\ 


Ho 


Ho (2> 

& : 


Reference 14 gives a tabulation of this function, which is rather 
important in the study of the aerodynamic phenomenon for the oscil- 
lating wing and which is reproduced in table 3. Kussner uses instead 
the function T correlated with \ by the relation 


T = 1 - 2 ,\ 


whereas in the paper by Kassner and Fingado (reference 15), the func- 
tion P = 1 - k is used with argument V = ui/2. 

American publications use the function C introduced by Theo- 
dorsen, which has the same definition as P. 

By means of equations (30) and (31), the pressures are easily 
computed for any type of oscillation (for example, translatory, rota- 
tional of the entire wing, or rotational of the flap). The coeffi- 
cients of the aerodynamic actions have been determined in various 
publications. A complete tabulation for the case of a wing with a 
flap hinged at the forward edge is found in reference 16, the com- 
putations for which were developed by the national institute for 
theoretical applications on the basis of the formulas of Kussner. 

In a recent publication by Kussner (reference 17), the case of the 
profile with a flap and with a tab hinged to the flap is treated. 

The fact that the hinges of the two movable parts can be retracted 
with respect to the corresponding leading edges is taken into account. 

16. Analysis of pressures on airfoil in motion in nonperturbed 
alr * - Equations (27) and (28) are considered, with the assumption 
that the values of w are due only to the motion of the wing and can 
therefore be expressed by equation (1). It is first assumed that V 
is constant. Then 


H 


w + 


v St 


dx 





dz 

St 


+ 



d 2 z 

St 2 


dx + 
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Inasmuch as adding a constant to the value of H does not 
change the result of equation ( 26 ) > 

x |ax-« 

and therefore 

H " v if + 2 sl + \ 

The first term in the second member together with the corres- 
ponding term contained in A q gives rise to the pressures that are 
denoted by Pq: 



n* 


itp 0 /2PV 2 ■ cot | 


dz 

s 


d£’ 




- sin & 



dz dd| 

5x cos - cos £ 


(33) 


where dz/ctoc is computed at the point of integration x' .* - cos . 
These pressures are those that are obtained if the wing in the actual 
configuration is under the conditions of steady motion.® 


The second term in the second member of equation (32) with the 
corresponding term contained in Aq. gives the pressures that are 
denoted by p^: 


jrpj/2pv = cot g 


0 






:dz d£* 
3t cos - 


cos ^ 


(34) 


These actions , which are proportional to the vertical velocities 
of the points of the wing, have the characteristic of damping forces 
and can, in part, be interpreted by kinematic considerations. Thus, 
if the wing is displaced without rotation with vertical velocity v, 

6 

In fact, substituting in this expression the value of u/V given 
by equation (37) of reference 11 for oz/dx 'and developing the computa- 
tion yields the value of "y = p/pv expressed by equation (31) pre- 
sented herein. 
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the intuitive result that the wing is in the same condition as if it 
ware at an angle of attack v/v is obtained from equation (34). 

Also, in regard to the effect of a torsional motion, a qualitative 
interpretation of can be given. If the wing is, for example, 

in diving rotation, it behaves with respect to the fluid as though 
it were curved upward. It is seen from such considerations that the 
focal moment that arises from this effect is turned in the opposite 
sense to the angular velocity (damping action). This consideration 
of the dynamic curvature has been put at the basis of the numerous 
approximate investigations on the aerodynamic coefficients of the 
vibrating wing (reference 9). These considerations would lead in 
substance to the computation of the values of p^ with the same rela- 
tion (equation (33)) that holds for p Q , in which dz/v dt is substi- 
tuted for dz/dx. This procedure leads to results that are quantita- 
tively in error because it is necessary to halve the second term in 
the expression of 


The third term of equation (32) gives rise to the pressures pg : 


On 


p g/ 2P = sin £ 


d$* 


COB d - COS 


0 


Ox' 


-1 


dt^ 


dx 


(35) 


These pressures, which result independently of the vertical 
accelerations of the points of the wing, represent an inertia effect 
of the mass of the circulating air. As evident from equation (35), 
the pressures pg do not depend directly on the local values of the 

accelerations, but on the entire distribution. Thus, in the case of 
the oscillation of a flap, although the forward part of the wing 
remains fixed, there are pressures over the entire chord. There can 
therefore, be no distribution of masses that are apparently capable 
of reproducing the inertia effect of the medium. 

The result of equation (35) is of interest for the rigid motion 
of the airfoil. In the case of translational oscillation, pressures 
are obtained that are distributed proportionately to sin d and give 
rise to the same resultant as though the mass of the cylinder of air 
circumscribed about the wing underwent the motion of the wing. It 
is easy in such a case to compute also the actions on a part of the 
chord. In the case of the rotational oscillation of the profile 
about the mean point, the pressures of inertia 3till correspond to a 
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system of masses distributed according to sin but the total mass 
in this case is equal to one-half that of the preceding case. The 
rigid wing thus undergoes in its motion an inertia action that can be 
thought of as reduced to two masses, each of a value equal to one- 
half of the mass of the circumscribed cylinder of air concentrated at 
a distance L/ V32" from either side of the mean point (always in 
regard to the computation of the resultant actions). 


It is observed that the pressures P2 have values that are 

independent of the velocity of advance and must therefore be sustained 
in air at rest. For such conditions, the results would not rigour - 
ously apply because the assumption of the smallness of the perturba- 
tions, as compared with V, does not hold. The results nevertheless 
agree with those that, for any particular case, have been obtained 
without the preceding assumptions and moreover apparently reproduce 
sufficiently well the actual phenomenon as it is found from some 
measurements by Cicala, in which the periods of the oscillation of 
a wing model in rarefied air and at normal pressure were compared. 

In the case of the phenomenon of the wing vibrations, the inertia 
pressures are not of great importance, whether because the additional 
masses represent a small part of the mass of the structure or because 
the previously described. actions, which exist independently of the 
velocity V, can be directly included in the computation if measure- 
ments of the mass of the structure are made by dynamic procedures. 

Equation (35), which permits computing the inertia pressures 
in closed form, is derived in reference 2, and the expression for p^ 

is contained in the same reference. The expression for pq presents 
a certain difference when compared with equation (34) in that it con- 
tains a term that in the preceding scheme is added to the transitory 
pressure. 

If the velocity V is not constant in the expression for H, 
there is derived from the term ^ P dx, in addition to the 

V J 5t ’ p x N 

quantity expressed by equation (32), the term dx, which 

because H is defined except for a constant is written z dv/dt. 
Hence, to the preceding computed pressures there are added the pres- 
sures pj, given by 


z d£* 
Icos £ - cos 


, dV 

P 3 / 2 P jjj = sin d 
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From this relation, there is obtained, for example, the following 
result: If the rectangular wing is displaced with velocity V not 

constant and with angle of attack a, which is constant, there is a 
force normal to the flight path represented hy the distribution 


dV 

P3 = 2pa dt Bln * 


which corresponds to the mass of the cylinder of air circumscribed 
about the wing subjected to the acceleration dv/dt. 

It is emphasized that the pressures f «* Po + Pl + ]?2 + i>3 are still 

added to the transitory pressures, depending on the values that the 
quantity w has assumed in the preceding instants. If the values 
of w are distributed linearly over the chord, w represents the 
value of w at the neutral rear point. The instantaneous pressures 
have the resultant passing through the focal point. 

The decomposition of the total pressures into instantaneous and 
transitory pressures has no absolute character in the sense that a 
part of the instantaneous pressures can be combined with the transi- 
tory (not vice versa, because a term containing the history of the 
motion is clearly distinguishable from the terms depending on the 
actual values). It nevertheless appears that the definition given 
herein is more natural because in the limiting case of steady motion 
the actions resulting from the group comprised of the instantaneous 
and the transitory pressures become zero. 

17. Remarks on treatment of unsteady motion of wing in two- 
dimensional case. - The first investigator to study the aerodynamic 
problem of the oscillating wing was Bimbaum, who made use of the 
concept indicated in Section 6 of the splitting of the circulation 
about the wing into bound and free components. Equation (7) is 
integrated for harmonic motion. By use of the complex variable 
notation, in this case. 


7 =7 (x) e*^ 

iQt 


e = c (x) e 
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As is easy to verify, equation (7) is obtained if 

e (x) - -we" u, ^J' X 7(x' ) e^*' dx' (36) 

•where 


u) * iQ/v 

Inasmuch as the free circulation is zero in correspondence with 
the leading edge of the profile (if no other sources of vortices 
occur upstream of the wing), in order that x = -1 and e = 0, the 
lower limit of the integral in equation (36) must he equal to -1. 
For the points of the wake, the integration is evidently limited to 
the chord, because 7 is zero outside the wing. 

On the basis of these results equation (8) assumes the form 



where x' and x” are variables of integration and w = WqS^* 1 . 

Bimbaum takes into account the condition of separation by 
expressing 7 by means of a combination of functions that satisfy 
this condition. The solution is sought in the form of a series expan- 
sion in powers of the reduced frequency QL/2V. The series converge 
rather slowly so that the results of Bimbaum are applicable only to 
rather low values of the reduced frequency. 

Wagner (reference 12) considered the problem of unsteady motion 
of the wing of infinite span. In the first part of reference 12, the 
treatment refers to the case in which w is constant over the entire 
chord but variable in time, a case of fundamental importance, as has 
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"been stated in Section 10. Making use of the condition of separation 
that imposes the circulation of the wing capable of giving rise to 
a finite velocity at the trailing edge, Wagner arrived at an integral 
equation that defines the circulation in the wake of the wing. With 
a procedure based on the moment theorems, Wagner gave an expression 
for the computation of the lift and of the moment on the wing on the 
basis of the circulation in the wake. In particular, the computa- 
tions were developed for the take-off motion of the wing. The case 
of rotation of the wing about a point of the chord was also considered. 


The problem of the unsteady motion of a wing was also considered 
by Glauert . In an initial paper (reference 18), he made use of the 
hypothesis that the circulation remains constant, a hypothesis that 
considerably limits the importance of the results. 7 In a succeed ing 
paper (reference 20), he took into account the variation of the 
circulation. In contrast to Bimbaum, Glauert sought to obtain 
directly the total circulation of the wing, which is divided into a 
part that would correspond to the case in which the vortex wake would 
be absent, and into a circulation induced by the vortices downstream 
of the wing. It is simple to show the equivalence of the relations 
assumed by Glauert for computing the lift P and the moment M with 
respect to the center point of the wing with the relations that are 
obtained on the basis of the bound circulation. From equation (4), 


p/pv = J"* (g-e) dx * K - J"* c 

m/pV = J (g-e ) x dx 


-\ 


dx 


> 

J 


(37) 


Considering that for x = -1, € * 0, and for x = 1, according 
to equation (3), Vc = - dK/dt, integration by parte yields 



—t 

Lamb (reference 19) also 
restriction. 


treated the problem with the same 
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With the aid of these relations, equations (37) become 


P/p 


w dK d 

+ dt + 3t 


*\ 


gx dx 


m/p = y 


« di + i « . i * 
^2 dt 2 3t 


-1 


n 

-i 




(38) 


dx 


considering that, according to equation (5), 


5x 5t 

and that the sign of the differentiation can he taken outside the 
integrals because the limits are independent of the time. Equa- 
tions (38), which are the equations used by Glauert, serve for the 
computation of the resulting aerodynamic actions but do not lend 
themselves directly to the determination of the pressure distribution, 
a computation that is necessary if it is desired to know the actions 
on part of the wing (flap). The equations obtained by Wagner are 
also subject to the same limitation. 

With the values of g expressed with the aid of the instantaneous 
characteristics of the rigid motion of the wing on the basis of the 
circulation existing in the wake, Glauert determined the lift and the 
moment with the aid of equations (38) and arrived at expressions 
agreeing with those of Wagners For harmonic motion, in which the 
circulation is distributed in the wake according to the sinusoidal 
law, the solution of the problem reduces to the determination of 
certain integrals that Glauert obtained by approximate procedures 
that limit the results to values not much higher than the reduced 
frequency. 

On the basis of the work of Bimbaum, Kussner (reference 8) again 
took up the problem of the oscillating wing, ass umi ng for the bound 
circulation the functions that are used for steady motion (refer- 
ence 11, p. 184). The corresponding w, expressed as a function of x, 
consists of a polynomial and of a trigonometric function multiplied 
by a factor depending on G 3, which Kussner computed by means of a 
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power series. After rather laborious computations, Klxssner obtained 
the pressure distribution on the wing corresponding to the motion of 
translational or rotational oscillation of the wing and also, approx- 
imately to the motion of the flap, for which he derived the coeffi- 
cients of the aerodynamic actions for the field of variation of G3, 
which is of interest for the phenomenon of wing vibration. 

Connected with the investigations of Wagner and Glauert is the 
theory developed by Theodorsen (reference 4), in which there are 
separately determined by use of the methods of conformal transfor- 
mation the potential function on the wing for simple fun dam ental 
motions of the wing with flap. From the potential function, Theodorsen 
determined the pressures with the aid of the Bernoulli’s equation 
generalized for nonsteady motion (reference 11, p. 38) : 


2 + l£ + = o 

P 2 3t 


(39) 


(where cp is the potential function). 


In the computation of the difference in the pressure between the 
lower and upper surfaces of the wing, it is necessary to consider, 
according to equation (39), the quantity derived from the variation 
with time of the difference in potential existing between the two sur- 
faces of the wing. This difference is measured from the circulation 
of the velocity, which is determined by following a path (fig. 3) that 
Joins the two points situated on the opposite surfaces and passes, 
always in the proximity of the wing, through the forward edge of the 


wing, thus obtaining the quantity 


II 


g dx; it is thus seen that 


the second term in the second member of equation (2) represents the 
corresponding term in dqp/dt in equation (39). 


Theodorsen* s treatment of the problem represents a narked advance 
with respect to the preceding work because it determines the distri- 
bution of the pressures on the wing and hence, in contrast to the work 
of Glauert and Wagner, permits the computation of all the coefficients 
of the aerodynamic actions for the wing with a flap. Also, because 
the integrals with which the effect of the vortical system of the wake 
is computed are solved by means of Hankel functions, all restrictions 
on the -value of the reduced frequency are thus eliminated by use of 
the existing tabulations*. 
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Independently of Theodorsen, and almost simultaneously. Cicala 
arrived at the solution of the problem of oscillatory motion of a 
deformable wing with any law by the same method followed by Bimbaum 
and Kiissner. In reference 9, it is proven that a class of functions 
exists, depending on the reduced frequency, that represents the dis- 
tribution of the bound vortices, which correspond to the velocity w 
distributed over the wing by a particularly simple law. The first 
of these functions, corresponding to constant w on the chord, con- 
tains the Hankel functions of the reduced frequency; the others are 
essentially represented by the y n of equation (14) and give rise to 

the values of w n in equation (12). In thiB manner, relations were 
obtained by means of which the coefficients of the Fourier series of 
the bound circulation on the wing can be computed in closed form as 
a function of the coefficients of the series for w. The coefficients 
of the aerodynamic actions for the wing with flap were thus computed. 

In a succeeding report (reference 2), it was also shown how the pres- 
sures depending on the second power of the reduced frequency (inertia 
pressures) and those proportional to the first power (pressures p^) 

could be computed in closed form without developing them into Fourier 
series. 

At the same time, Kussner arrived at the solution of reference 9 
by a procedure described in reference 21, some results of which were 
anticipated in reference 22. The general case of nonsteady motion 
was also treated in reference 21, where the discontinuous motion was 
studied on the basis of the solution for the harmonic motion with the 
use of the Fourier integral. The case of a stationary gust was studied 
and tests were conducted (reference 23) confirming the result obtained 
that the focal moment remains zero during passage through the gust. 

The solution of Kussner was used by Dietze (reference 24) to compute 
the resultant of the actions on the flap (in the preceding papers 
only the hinge moments were computed); it was also used by Krall 
(reference 16) to elaborate, with the aid of the National Institute 
for Applied Computations, the tables of the aerodynamic coefficients 
for the oscillating wing, and was used by Dietze again (reference 25) 
for the computation of the coefficients for the wing with a flap and 
a tab hinged to the flap. 

Kassner and Fingado (reference 15) also succeeded in computing 
the actions corresponding to the oscillatory rigid motion of the wing, 
making use of the Wagner's expressions by which, in the case of har- 
monic motion, the evaluation of the integrals relative to the effects 
of the vortices in the wake was investigated with the aid of Hankel 
functions by Borb&ly (reference 26). With the aid of the Wagner's 
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expressions, Ellenberger computed, the resultant actions on the wing for 
a flap rotating according to a general function of time (reference 27 ) . 

In a summarizing note, Jaeckel (reference 28) established a coor- 
dination between the procedures of Glauert, of Lamb, and of Birnbaum- 
Kiissner for the solution of unsteady motion of wings and considered 
also the case of the wing with a variable chord. Jaeckel also pub- 
lished a systematic derivation (reference 29) of the results that were 
given with rather synthetic justification by Kussner . A clear review 
of the theories on the wing in unsteady motion was given by Lyon 
(reference 30). 

The results of the preceding studies, with some. further develop- 
ment, are treated by von Kannan and Sears (reference 7); a derivation 
procedure is developed that presents in an intuitive form the mathemati- 
cal fundamentals of the investigation. The computation is restricted to 
the determination of the lift and the aerodynamic moment, which are 
computed with the aid of the expressions, respectively, 

p/p=^2rixi (40- ) 

M/P - 5 I Tj X\ (40”) 


The second member of the first expression represents the rate of 
variation of the moment of circulation of the system measured with 
respect to a fixed point; in an interval of time dt, this moment 
should vary as the bound vortices are displaced by the amount V dt, 
while the position of the free vortices and therefore the moment with 
respect to the fixed axis have not changed. The total variation will 
therefore be given by the product of the total bound circulation and 
V dt and therefore, when multiplied by the density and divided by dt, 
will give the lift. By analogous reasoning, the second of equa- 
tions (40) is verified. Thus, the square of the distance from the 
fixed point varies by the amount 2 VI dt for the bound vortices while 
it remains constant for the free vortices. Hence, the variations of 
the second member of equation (40" ) will be represented by the moment 
of the bound circulation or the aerodynamic moment except for the 
factor p . 

From these relations, von Ka'rman and Sears derived expressions for 
the lift and the aerodynamic moment, which present a better generali- 
zation than those of Glauert inaemuch as they can also be applied to 



42 


MCA TM 1277 


the deformable profile; but the relations still do not permit compu- 
tations on parts of the wing. The authors divided the total circula- 
tion into a part that would he obtained in the absence of the vortex 
wake (so-called quasi- stationary system) and into the induction of the 
vortices of the wake. The case of the general motion of the profile 
and of a stationary gust is also considered. 

Garrick (references 31 and 32) brought out the relation (Laplace 
transform) that exists between the Wagner function, which gives the 
circulation by the elementary discontinuity, and the function that 
holds for the harmonic motion and proposed the use of approximate 
expressions to represent the function of Wagner, which would be use- 
ful in the analytic solution of various problems of unsteady motion. 

Possio, making use of the acceleration potential, considered the 
problem of the discontinuous motion of a wing (reference 10) and the 
case of the stationary gust. 

In a recent publication (reference 17), Kussner and Schwarz indi- 
cate the relations with which the pressures on the profile, can be 
computed without making use of Fourier series but using integrations. 
These equations were applied to the determination of the aerodynamic 
coefficients for the wing with flap and hinged tab on the flap, when 
the case is considered in which the hinges are set back relative to 
the leading edge of the moving parts. The same relation for the com- 
putation of the pressures on the oscillating wing was derived in a 
different manner in a report by Schwarz (reference 5), in which a 
clear and rigorous derivation of the known solution of the aerodynamic 
problem of the oscillating wing is developed. 

In reference 13, Kussner gives the general solution of unsteady 
motion in the two-dimensional case and the functions of Wagner end 
those relative to the aerodynamic actions produced by a gust are com- 
puted with greater precision. 

The general case of unsteady motion has also been treated in a 
report by Sohngen (reference 6). The solution is put in substantially 
the form indicated in Section 13. 
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PART II. - WING OF FINITE ASPECT RATIO 

The relations that connect the velocities induced hy the vor- 
ticity of the inductor system of a wing of finite aspect ratio are 
derived herein and a form for these expressions is sought that lends 
itself to a future refinement of the investigation that, up to the 
present, has "been conducted with approximate procedures. An evalua- 
tion is then given of the various approximations that have "been used 
in the theory of the finite wing in unsteady motion. 

The consideration is restricted to essentially rectilinear 
motion of the lifting surface; that is, it is assumed that the velo- 
city of the points of the wing give small deviations relative to a 
mean value V, which maintains its direction unchanged with respect 
to the motionless fluid at infinity. It is assumed that V is small 
compared with the velocity of sound. The origin of the orthogonal 
axes x, y, and z is located at the point at which the induced velo- 
city is measured (point of induction). The x-axis is taken parallel 
to V and in opposite direction; the y-axis is normal to the x-axis 
and is contained in the fundamental plane that, during the motion, is 
at a very small distance from the points of the lifting wing, which 
is assumed to "be of infinitesimal thickness and curvature; the z -axis 
is perpendicular to the x- and y-axes and directed downwards. The 
vortices having an axis parallel to the x-axis (longitudinal vortices) 
and to the y-axis (transverse vortices) are considered positive if 
turned in the sense that carries the positive directions of x and y 
on z. In the middle section of the wing is located the origin of the 
£- and t) - axes parallel to x and y. The semispan of the wing is 
denoted by b, so that for the points -b < q < b, i 0 and t) 0 denote 
the coordinates of the point of induction, hence 


x 


e - e 


0 


y - n - no 


Let £n(ri ) *e the equation of the leading edge; £ p (q) the equa- 
tion of the trailing edge; i v 4 n L the chord; ft the pulsation; 
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u ) the imaginary factor iQ/Vj u> the reduced frequency Ql/2V; 

the integral taken over the chord from the leading to the trailing 
edge. The other symbols are defined in the text or in Part I. 

18. Inductor elements in tridimensional case. - Let P (fig. 8) 

"be a point of the wing that moves relative to the fluid with vector 
velocity y of magnitude v and having any direction; the linear ele- 
ment dl through P normal to y supports the "bound vorticity T; 
that is, supports the aerodynamic action pvr dl. In the interval 
dt that precedes the actual instant, the point P starting from P' 
is displaced by the segment da = v dt. The total inductor is then 
changed in that the bound vorticity element has come to occupy the 
position 1-2 from the position 3-4, leaving behind it free vortices 
and simultaneously creating the two longitudinal elements 2-3 and 1-4 
of equal circulation T . There is thus added to the preexisting 
inductor system a closed vortex element 1-2-3-4; on the side 3-4, which 
at the instant t there exists the vorticity element of intensity 
-dT liberated from the bound vortex, it may be assumed that there 

clX " 1 

simultaneously exist the vorticity r - — do, which existed at the 

time t-dt and the element of intensity T with oppositely directed 
sign constituting the fourth side of the circuit 1-2-3-4. 

With the aid of the formula of Biot, it is found by simple com- 
putation that the closed vortex element induces, at a point 0 at 
distance r from P, a velocity that, except for infinitesimals of 
the higher order, may be written, as® 

d^w = - T dl dj/4«r^ 

Whatever happened in the interval dt occurred in all the pre- 
ceding time starting from the instant at which a force on the lifting 
element has arisen; that is, since an element of the bound vorticity 
was created. In the problems that ordinarily present themselves, it 
is assumed that v always has the same direction parallel to the 
x-axis, or more precisely, the deviations with respect to this direc- 
tion are considered sufficiently small. Let do and dl therefore 

O p 

The symbol d is used to indicate a quantity that, divided by 
the product A l AO, has a finite limit when these elements approach zero. 
With the assumed signs, w is directed upward and is therefore 
negative . 
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"be parallel to x and y, respectively. The symbol o denotes the 
distance measured from the point P of coordinates z and y, 
parallel to x, assuming that the circulation T on the element of 
the lifting surface considered is expressed as a function of Q* 
Summing the effects of the closed vortex elements created in the 
preceding time, the velocity at 0 is obtained: 


dv 


dy r do 
4* I r 3 


d£ 

4k 


Tdo 


/\/(x+o) 2 + y2 3 


(41) 


The velocity induced by the bound-vortex element and by the 
system generated by it for r ^ 0 is identical to that which is 
obtained on the basis of the concept of the acceleration poten- 
tial and which is defined as the field of the pressure dipoles 
(reference 33). 

The particular case of steady motion ( T constant) is con- 
sidered. Inasmuch as 



a quantity that is denoted by f^c+o^y), then from equation (41), 

- 4k dw « Tf(x,y) dy (43) 

where dw denotes the velocity induced by the elementary horseshoe 
vortex of frontal Bide dy and circulation T . The velocity (fig. 9b) 

induced by the semi vortex 2 «» is i (1 - x 

7 ' V * 2 + y 2 

summing the velocity corresponding to the other semivortex, which 
has the opposite sign to the first, there is obtained, except for 
infinitesimals of higher order, the quantity 


) 


and therefore. 
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Adding to this quantity the velocity corresponding to the 
frontal sequent yields the value of dw given in equation (43). 
The equivalence, which in this case holds between the inductor 
system of the pressure dipole and the elementary horseshoe vortex, 
can also be seen by considering that the sides parallel to y of 
the continuous vortex circuits (indicated by the small circles in 
fig. 9(a)) are canceled if r has the same value for all. This 
inductor element is called a bivortex and a section of it closed 
in front and behind is called a segment of a bi vortex. 

Making vise of equation (42) and integrating by parts yields, 
from equation (41) , 


- 4n dw 


r f (x,y) dy + dy 


dT 

do 


f(x+o,y) d 


(44) 


The operation that leads from equation (41) to equation (44) 
transforms the inductor system constituted by the closed circuits 
into a system of bivortices: a generating bivortex originating 

on the wing, denoted as a bound bi vortex, and a row of free vortices. 
In a distance do of the wake of the bound bivortex, free bivorticeB 
originate the circulation of which is represented by the variation 
that the intensity of the bound bivortex has undergone in passing 

through the distance do; that is, ^ d(7. 

do 


For harmonic motion, • 


r - Fe iQt 


where F is constant and 0 is the frequency of the motion. In 
this case, V is always assumed constant. The circulation of 
element 12 at the time in which the point P of the wing was set 
back with respect to the actual position by the amount is repre- 
sented by 

T(o) - 

that is, with respect to the actual circulation I" 1 , a lagging phase 
shift represented by iftj/v = do. Making use of this result and 
denoting by #r dy the velocity induced by the bound bivortex and 
by the system of the wake in the case of harmonic motion yield 
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^as 


- 4nw « 


e-^do 


VCx-KJ) 2 + y 2 


(45) 


By setting 


x * Q x/v y = ft y/v o - Qo/v 


equation (45) can "be written in the form 


4jcw a (i? $ 


(46) 


where $ is the function of two variables x and y defined by the 
relation * 


$ 



(46a) 


For x = 0, this function has been computed by Muller (refer- 
ence 34, p» 36). For x ^ Of 


$ (x,y) « e lx 

®(0 ,y) - 

o? n 

-iu . 

© du 


- C 

0 (u 2 +y 2 ) 3 / 2 


With the aid of this relation, which is easily obtained from 
equation (46a), the computation is reduced to the tabulated function 
and to an integration between finite limits that can be carried out 

by graphical or numerical methods. The computation of was carried 
out for y » 2. When u 


$ * $>» + i$ n 

there are plotted as abscissas and ordinates in figure 10 the real 
and imaginary parts $• and of the function and the curve with 
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the values of x drawn. This curve serves for the computation of 
the velocity induced "by the pressure dipole at the points' of a sec- 
tion at distance y * 2v/^ from the inductor elements. The point 
x = 0 is at the origin of the dipole. The points that are found 
upstream correspond to the positive values of x. For this segment, 
the values of w decrease continuously hut all have slightly dif- 
ferent phases. Proceeding toward the downstream region, the points 
shifted hack with respect to the front of the hivortex generator 
(x < 0) encounter velocities that vary little in magnitude hut rela- 
tively more in phase. It is evident that at a great distance from 
the origin the velocities vary on the parallel considered hy a sinu- 
soidal law and therefore the representative point of figure 10 for 
negative values of x increasing in absolute value will approach a 
certain limiting circle with center at the origin. 

19. Correction of divergent expressions. - The expressions of 
the preceding section are sufficiently well adopted for the computa- 
tion of the velocities corresponding to the inductor system of a wing. 
Special methods are required, however, for the pointB of infinity that 
are presented hy the functions under the integral sign present. 

The case of steady motion is first considered. A hound vortex 
filament AB (fig. 11(a)) of circulation V varying from point to 
point is considered and with this filament is associated the system 
of longitudinal vortices that is ehed from the points of AB according 
to the condition of Helmholtz. The velocity that these elements 
induced at a point 0, which is outside the vortex wake of the fila- 
ment AB, can he computed hy means of the relation 

PA 

- 4jtw * T f (x,y) dy (47) 

JB 

where x and y are the coordinates of a point of the filament AB 
at which T is the circulation. This expression reduces the inductor 
elements to a system of hivortices that originate on the filament AB 
(scheme of fig. 11(h) for A— > 0). The demonstration of this equa- 
tion is given in reference 35. An intuitive Justification is obtained 
from figure 11 (b). The segments of the dotted vortices in the limit 
reproduce the effect of the generator filament AB, whereas the other 
longitudinal elements constitute the system of the marginal vortices; 
each pair of contigous elements function as a single inductor ele- 
ment of circu la tion equal to the difference of the circulations of 
the two elements that compose it. 
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This expression is applied to an indefinite vortex of unit cir- 
culation at a distance x Q from the point of induction. If x Q is 
positive, f always remains finite. In fact. 



If iq is negative, f increases indefinitely for y— K). 

Hence, in a singular integral, it is necessary in computing equa- 
tion (47) to exclude from the integration the s m all segment from 
y = - 5 to y « +5 and then make 5 approach zero. The expression 
that is thus obtained, however, has no limit, as is intuitively evi- 
dent from an examination of figure 12. The bivortices (fig* 11(a)), 
into which equation (47) transforms the inductor system when the ele- 
ments contained in strip 26 (which Includes the point 0) are 
excluded, are equivalent to two angular vortices (fig. 12(b)) (because 
the semi vortices Joined by the small circles are canceled) and these 
elements, on approaching 0, induce a velocity that increases without 
limit. For this reason, the integral equation (47) diverges if the 
abscissa xq from the point at which the filament AB cuts the x-axis 

is negative, as in the simplified case of figure 12; in the computa- 
tion of the principal value of the singular integral, an expression 
is obtained that increases indefinitely and that represents the velo- 
city induced by the two semi vortices, which are indefinitely removed 
from the point 0. (See fig. 12.) 

In order to eliminate the previously discussed divergence, the 
following artifice may be applied: There are added to the elements 

of the filament AB those of the indefinite vortex I (fig. 11) 
passing through the point Q at which the x-axis through the point 0 
intersects the filament AB and having the circulation -Po, where 
r 0 is the value of T at Q. The system thus obtained gives a velo- 
city at 0 that is denoted by Dw and, according to equation (47), 
is represented by 

4« dw * jr 0 f(xo,y) - (48) 

U- 00 


The first term of this integral represents the effects of the 
bivortices having their origin on vortex I; the second term refers 
to the bivortices having origins on filament AB (hence the integra- 
tion is extended from y of the point A to that of B). The two 
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terms are divergent when taken separately, hut give a correct expres- 
sion when added. The intuitive reasoning for this lies in the fact 
that the two semivortices that, according to the scheme of fig- 
ure 12, include the point 0 and arise from the vortex I find 
compensating terms in the elements arising from the filament AB. 

The velocity induced hy 'the vortex I is then added to Dw. 

In the case of a system 2 of inductor elements analogous to 
the filament AB with corresponding marginal vortices, the following 
procedure may he used: 

I) The transverse vortices cut out from the section hy the 
induced points are prolonged indefinitely, thus obtaining a system 

the induced velocity of which w* is computed hy the formulas 

of the two-dimensional motion. 

II) The quantity Dw, which constitutes the velocity induced 
hy the system 1 2 * 2 - Ii obtained hy superposing on the inductor 
system Z the system Zi with reversed sign, is computed with the 
aid of equation (47) .and added to w* . 

The advantage of the preceding procedure lies in the fact that 
both w* and Dw have a homogeneous inductor system (indefinite 
vortices for 2^, 'bivortices for Zg) an ^ therefore free or hound 

longitudinal and transverse vortices need not he separately considered. 
The advantage is reflected in the simplicity of the formulation of w 
and in the rapidity with which the results of the existing theories 
for the approximate solution of the problem under consideration are 
derived from this scheme. 

The same procedure can he applied if the motion is unsteady. 

The hound -vortex system impresses on the fluid an imprint represented 
hy a similar vortex configuration corresponding to the variation of 
the intensity of the generating system. These vortex systems carried 
hy the stream may he treated as shown for the hound system. The 
plane system 2^ including the free vortices can he analyzed hy the 

formulae of part I. The system I>z is made up of hound vortices 
with the corresponding wake of the row of free elements (that is, 
pressure dipoles) and can therefore he analyzed with the aid of 
equation (48). 

20. Vortex system of wing. - The velocities induced hy the vortex 
system of a wing can thus he computed hy adding to w' , induced hy 
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the plane system that would obtain if the transverse vortices cut 
out from the section through the point of induction were indefinitely 
extended, the velocities Dw that, as a function of the hound vor- 
tices, are written as 


Dw => l tv di dy - 7 n $ dx dy ( 49 ) 

Js Js 1 

If harmonic motion is considered, the functions w can he 

derived from equations (45) and (46). 9 If the intensity of the gen- 
erator bivortex varies according to a different law, it is necessary 
to make use of equation (41), for which the function T(a) must 
he known. From this relation the function b « dw/dy T(0) is 
obtained, which, in general, varies with time. 

The first of the two integrals of equation (49) is extended to 
the surface S of the wing, y being the bound vorticity, the 
second integral is extended to the strip included between the 

lines Z]_ and t]_ parallel to y passing through the points at 

which the chord through the point 0 cuts the leading edge l and 
trailing edge t; the vorticity y Q on this strip is that which is 

on the wing at the section through 0. On each element of the sur- 
faces S and there originates a bound bivortex of circulation 

7 dx dy connected with the pressure, which acts on the element by 
the relation 


p ss PV7 

In the computation of Dw, as in the computation of w* , the 
total circulation g along y can be referred to instead of 7 , 
the two quantities being connected by the relations of part I. When 
the effects of the elements of a section y « constant are considered, 
a system of bivortices will be obtained the origins of which are dis- 
tributed either on the chord or in the wake, the intensity of the dis- 
tribution being represented by g. The expressions are restricted to 
harmonic motion and are written as 


9 A 

As the reduced frequency uj increases, the modulus of w 

decreases and therefore also Dw, which represent the correction 

velocities due to the finite span. 
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g ■ gC^n) 6 
K - gdx - K(ti) e 1 ^ 

In the -wake at a distance 0 from the rear edge in a segment do, 
bivortices originate that have the circulation 


if^Ldo 

do 

equ al to the change that the total circulation K has undergone in 
the interval in which the rear edge has passed through the distance do 
at the time t-o/v in which K(o) was displaced with respect to the 
actual value "by Qo /v : 


K( 0 ) « K e" wa 


The velocity induced hy the elements of the section considered 
is therefore proportional to 


r,«.r 

Jl Jo 


gf dx + 1 d ( Ke ~ axr ) . f do - r gf dx -ojk r 

1 d ° Jl Jo 


e - u)Of da ( 50 ) 


The symbol ip denotes the abscissa of the rear edge of the 
chord (that is, 2p ■ Ip - £ q)> therefore 


f ' ' (x,y) « t (x,y) - f (x P ,y) 

It is easily found by using equations (44), that equation (46) 
can be written in the form 


- u, 2 4 ■ f(x,y) 


u> I f(x+o,y) e“ w0 do 

Jo 
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Hence, equation (50) can "be transformed to 


l"' 


dx - U) Z K$ (x p ,y) 


(50* ) 


The expression for the velocities Dw at the point £ q, tj q is 
therefore obtained : 


- 4tr Dw = 


r\r 

J -b J L(t] ) 


gU>n) d £ 


dTJ . s(bno) f ' (x,y) d£ 


IL(rio) 


0) 


! j b $(xp,y) * K(t]) dri +0) 2 I ^(xo,y) • K(t) 0 ) d^ 

J •!) J -OD 


(51) 


indicates the integral taken from |„(ti) and 
the quantity Q p (ti 0 ) - Q/v. 

The first of the integrals in the second member of equation (51) 
represents the induction of a system of segments of hivortices having 
their origin on the points of the wing and ending on the rear edge of 
the chord, and which are of intensity equal to the total vorticity at 
the point of origin. The second integral represents the effects of 
the analogous inductor system existing in the strip S^. The third 

integral expresses the induction of the wake system constituted by a 
system of pressure dipoles with origin on the trailing edge of the 
wing. The fourth integral refers to the elements of the wake having 
their origin on the line t 1 « 

If 


Where X(Tl) 

and with Xq 
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the first two terms in the second member of equation (51), which 
corresponds to the velocity denoted by w 2 , can be transformed by 

integration by parts so as to assume the form 


4jtwo = 


pb 

* 1 

g( ) dc, 

^00 

dri 

y 2 

-b J 

r 3 

L(ri) J 

y 2 

-« J 


G(|>Tq) dT) 


L(tl 0 ) 


where 


r 


=/^x 2 +y 2 


It is then observed that, in regard to the effects of the ele- 
ments of the wake, the difficulties arising from the divergence of 
eqmtion (47 ) cannot appear because all the points of the wing are 
located upstream of the line from which the system of the wake 
originates. Hence, for these elements the decomposition into the 
systems and tz can be avoided and the total w = w* + Dw can 

be directly computed. The velocity denoted by W 3 is thus obtained 
and is expressed by 


4 ™3 


Pb 

2 $(xp>y) 

J-b 


K(t]) dr) 


2 K(n 0 ) 


If W]_ denotes the velocities induced by the plane system that 
is obtained by indefinitely prolonging the total vorticities cut by 
the section containing the point 0 ; that is, if 


2jrw^ 


P g(5>io> 


Win) 


X 


the resultant velocities are expressed by the sum w « w^ + w g + w^. 
In general, w^ must represent the preponderant part of w, which 

facilitates development of the methods of iteration for the computa- 
tion of the function g for an assigned w. The velocities w 2 , 

which represent the part whose computation presents the greatest 
analytical complexities, are independent of the reduced frequency and 
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are therefore determined, only once for any frequency of oscillation. 
The values of W3, which depend on the frequency of the motion, are 

expressed hy simple integrals and are, moreover, zero for all the 
distributions of the circulation (equation (11)) for which the 
integral K is zero. 

21. Approximate theories. - A first approximate treatment of 
the problem of unsteady motion of the finite wing was developed by 
Cicala (reference 36). The principle of approximation there assumed 
finds simple formulation if reference is made to the expressions of 
the preceding section. Although use is made of the exact solution 
regarding w' , in the computation of Dw the segments of the bi vor- 
tices of to which the first two terms of equation (51) correspond, 

are neglected and, moreover, it is assumed that on every chord the 
velocity induced by the system of the wake is constant and equal to 
that which would obtain if the system started at the point of induc- 
tion. Hence, referring to equation (51), in addition to neglecting 
the first two terms of the second member, in the computation of the 
other two by $, a function of the abscissa and of the ordinate of 
the point of t with respect to 0, the value corresponding to 
x = 0 is substituted. As a consequence of the first approximation, 

Dw = 0 if K = 0. Hence, if g is represented by a linear combi- 
nation of gn of equation (11), the corresponding velocities are the 
w n of equation (12) and the pressures are obtained from the y n 
defined by equation (14). Reciprocally, if the vertical velocities 
on the various chords can be represented by a combination of w n , 

•the corresponding circulations and pressures can be computed on the 
various chords, section by section, as though the motion were two 
dimensional. This fundamental simplification permits reducing the 
tridimensional problem to a single case; for example, that of ver- 
bical velocities constant on the different chords: any distribution 
of w that is developed in the series of equation (16) (in this 
case, A n is a function of the coordinate measured along the 
span) requires particular examination only for the circulation 
corresponding to the w defined by equation (17) (this simplifica- 
tion is also, in general, variable from one section to the next). 

_The bound circulation corresponding to the velocity 
w = w = constant over the entire chord is obtained on the basis of 
equation (19) in which, if harmonic motion is considered, according 
to equation (22) and the results of Section 15, the following rela- 
tion must be substituted: 


C = w(l-\) 
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From 7 , e is then computed with the aid of equation (36) 
and then the total circulation yields 


K « J"* (y+e) dx 

VThen the computations are made, it is found that "between w 
and K or "between the total amplitudes of the quantities W and K 
there exists the relation 

w (H^ 0 “ ei w K - 2WL 


where H is the Bankel function of parameter w . 

The velocity W is that induced "by the system to this 
quantity is added the velocity due to the system 2 2 , which, for the 

assumed simplifications, is also constant on every chord. With the 
computations developed in reference 35, equation (51) is transformed 
into the relation 


- 4 * V" - J_ b if N dT l (53) 


where 


"H coordinate measured from middle section of wing parallel 
to y (which is measured from point 0 ) 


N = — - (uF 

y 


F function of variable y « fiy/v defined "by 
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This function is computed by this relation for y > 0. For 
y < o, F(y) = -F(-y). The f motion F is tabulated in reference 37; 
in a recent report by Kussner (reference 33)* the function N is 
tabulated* 

Adding the -velocities induced by the systems £ i and yields 


_ u)(H ( 2) 0 )- iH^i) e^K 

W " 2L 



(54) 


This equation* on the basis of known -values of W* defiiies the 
distribution of the total circulations along the span. For u) ■ 0, 
this reduces to the lntegrodlfferential equation of Prandtl* Inasmuch 
as* for the simplifications assumed* the distribution of 7 over the 
chords is similar to that of the two-dimensional motion* on the basis 
of the values of the total circulations* the problem is completely 
solved. Given the series of equation (16) on w* the circulations 
and the pressures corresponding to the part that is developable in 
the series of vn is computed as if the motion were two-dimensional* 
whereas for the remainder w* the circulations are distributed on the 
chord as if the motion were two-dimensional and along the span of the 
basis of the solution of equation (54). 


In a succeeding note (reference 37), the procedure was applied 
to the determination of the aerodynamic coefficients for the oscil- 
lating wing. ' 

Independently of reference 36, Borb&ly proposed a type of approxi- 
mation for the computation of w for the finite wing (reference 38). 
Reference 35 shows that the expressions that-Borb&ly elaborated for 
the computation of Dw in the particular case of the elliptical 
distribution of K along the span agree with the results that* for 
this case, were derived on the basis of equation (53). 


Possio* concerned with the problem of the stability of small 
oscillations of the wing considered as a rigid body, also analyzed 
the problem of the oscillating wing of finite span. Making use of 
the concept of acceleration potential* he derived, equations (9) and 
(10), which* however, as has been shown in Section 8* are also Justi- 
fiable on the basis of the concepts of vortices. The solution is 
expressed by Possio in the form of a power series of the parameter 
Qb/v . The value limited by this parameter and the smallness of the 
ratio L/b (large aspect ratios) in the series containing the powers 
of QL/2V Justify, for the computation of the w* corresponding to 
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the system Iq, the Introduction of simplifications that were not 

adopted in references 36 and 37. The assumed simplification in the 
computation of Dw, when it is reduced to the scheme represented hy 
equation (49), can he thus defined: The function "w is computed 

hy means of equation (46), substituting for $ the value that, 
according to equation (46a), corresponds to x * 0 and that is there- 
fore constant for the elements of each chord. Reference 39 contains 
the principles of the procedure. Some of the results are described 
in reference 40 and in greater detail in reference 41. More general 
cases of the mbtion are considered in references 42 to 44. Refer- 
ence 44 analyses the law of variation of the lif-fe on a rigid wing of 
elliptical plan form during the start of the motion (the same problem 
that was studied by Wagner in the two-dimensional field). In the 
computation of the velocities induced by the transverse vortices, 
there was assumed (in the simplification of equation (14) of refer- 
ence 44) an approximation different from that used in reference 43; 
the approximation gives for the case considered a greater precision 
of the results. 

Sears (reference 45) also studied the problem of the oscillating 
wing of rectangular plan form with approximate procedures. Reducing 
the computations to the scheme of the preceding section, the simpli- 
fication assumed consists of the suppression of W 2 , while a rigorous 
computation is proposed for w 3 , and mating use of the results tha t 

are obtained for the infinite wing with sinusoidal distribution of the 
circulations along the span. The computation, which is intended to 
eli mina te the errors inherent in the approximate theory of refer- 
ences 36 and 37 criticized by Sears, does not give results more accu- 
rate th a n that theory. As an example ,_the effects of the inductor 
elements of a chord L, for the case » 1 10 , are considered with 
the point of induction at a distance y * L. Apart from factors that 
need not be considered in a comparison, the induction can be expressed 
by the equation (50* ), which can be written as 



where $ ^ is the nondimensional factor 


10 


In the cases of wing vibrations encountered in practice, the 
reduced frequency turns about this value. 
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•(x,y) -*'! + »" x 

x "being the abscissa of the point, with the vorticity g and xp 
the abscissa of the rear end of the chord measured with respect to 
the position of the point of induction. The quantity $1 was com- 
puted for the elements of the circulation g at the leading edge, 
at the middle, and at the rear edge of the chord considered, and in 
the three cases, on modification of the position of the induction 
point, there were obtained (see fig. 13, in which are drawn as the 
abscissas and ordinates the real and imaginary parts of $ 1 , 
respectively) the curves I, II, and III in the figure with the values 
of I * xp/L. If the approximation that was made in reference 36 
is assumed, in the three cases, for any position of the induction 
point, the end of the representative vector is the point indicated 
by the double circle; with the aid of the approximation of Sears 
for all three cases, the representative points are those of curve III. 
The error with either approximation is large. If the greater com- 
plexity of the computations required by the solution of Sears is 
considered, the advantage of his approximation is questionable. 

In a recent publication (reference 33), Kussner, wwictng use 
of the acceleration potential, developed a new approximate theory 
of the oscillating wing of finite span. The approximation assumed 
is easily related to the expressions of the preceding section; 
referring to equation (49), the value of w given by equation (45), 
setting u * x + o, can be written in the form 


^U^y) = y 2 



(45*) 


Kussner* s solution can be obtained by setting the lower limit 
of this integral equal to zero. The Dw, on the "basis of this 

assumption, are proportional to e” w £ (where £ is the coordinate 
measured parallel to x from a fixed origin of the wing in an arbi- 
trary position); that is, the values of Dw are distributed by the 
sinusoidal law over each chord. If the distribution of the vorticity 
is considered, which on the basis of the solution of the two-dimensional 
motion corresponds to the velocities distributed according to thi s 
law, and if for this case the values of w of the system are 
summed, velocities are obtained that can be represented by the 
expression 
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rr iQt-w£ t -r -U>(£-Vt) 
w = We = We ' ' 


(55) 


where W is a function of q. 


By obtaining, on the basis of the solution of the two-dimensional 
problem, the velocities induced by the system 2q and adding the 
values of Dv computed with the aid of equation (49), simplified 
according to the preceding assumption with respect to the limit of 
the integral in equation (45) and transformed by operations analogous 
to those that led to equation (54) (see reference 35), the final equa- 
tion is obtained 


u>(h ( 2 ) 0 - iH^^jK* i 
" “ 2L(j 0 -iJi) 4* 

I 

in which H and J are cylinder functions of the parameter SJ and. 



r d 


d£* 
-b dT l 


N dq 


(56) 


This quantity has the same modulus as the total circulation K, 
but has a certain phase displacement with respect to it. 

By solving the integrodifferential equation (56), which has the 
same kernel as equation (54), the law of distribution of the vorticity 
on the wing corresponding to the •values of w given by equation (55) 
can be obtained. If it is assumed that i -Vt is the abscissa 
measured in a system of reference fixed with respect to the fluid, it 
is concluded that equation (55) represents a distribution of velocity 
having local values constant in time, as would be the case of dis- 
turbed air that presents, along the trajectory of the wing, vertical 
currents of constant velocities in time (stationary gust). The solu- 
tion indicated refers to the case of the stationary gust of sinusoidal 
form. The more general case of motion can be studied with the aid of 
the preceding solution and the solution of the two-dimensional problem 
when it is considered that, according to the approximation of Kiissner 
(as with the procedure of references 36 and 37), the values of Dw 
are zero for the distribution of the g that gives rise to zero 
values of the integral K and therefore of E' , which has the same 
modulus as E. Hence, in this case, to the values of w represented 
by a combination of w n expressed by equation (12) there correspond 
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the -values of g expressed by equation (11), as if the motion were 
two-dimensional. On this basis, reference 35 indicates the extension 
to general cases of the solution based on equation (56). Kussner 
makes the generalization by a different principle, which leads, 
however, to expressions that, in the limiting case of the infinite 
wing, give the exact solution already known. 

The approximate theory of Kussner therefore cannot be derived on • 
the basis of the vortices concept developed in Section 18, as is true 
for the solution proposed by Cicala, which is criticized in the pre- 
ceding note by Kussner as presenting arbitrary assumptions; the disa- 
greement between equations (54) and (56) finds its Justification in 
the different principle of approximation rather than in a fault of the 
derivation method based on the vortex concepts, which, according to 
Kussner, would lead to erroneous results. The approximations thus far 
assumed all lead to a somewhat inexact value of the induced velocities, 
as is shown in reference 35; the theories are all, except that of 
Sears, constructed so as to converge in the case of steady motion to 
the theory of the vortex filament of Prandtl, whose approximation has 
thus far been proven sufficient. On analytically examining some local 
values of the errors committed, ^ all are shown to be of little pre- 
cision, from the simplest to that of Sears, which consists of the most 
laborious application, or that of Kussner, which is based on the 
concept of pressure dipoles. Only the fact that in a limiting case 


11 If the solution is expressed in the form of a decreasing power 
series of the aspect ratio \, it is found that the error in the 
existing theories starts from the term in logX /X 2 . In reference 35, 
rather than analyzing the order of magnitude of the error, it is pre- 
ferred to carry out the computation for concrete cases so as to be 
able to compare the various approximations. 

The comparison is particularly evident by making reference to the 
concept of pressure dipoles. According to the principle followed in 
reference 39, the dipoles of the system £ 2 are transported parallel 
to the direction of the x-axis up to the induction point; according 
to the principle adopted in reference 33, these dipoles are given the 
same displacment and, in addition, the same phase displacment 
Q x/V (x is the abscissa of the dipole with respect to the induc- 
tion point); according to references 36 and 37, these elements are 
given the same displacement and phase shift Q(x-xp)/v. It is shown 
in reference 35 that the three approximations alter somewhat the 
value of the velocity produced by the pressure dipole. The affinity 
of the three principles is evident. 
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the approximations lead to a theory that has shown itself satis- 
factory in applications indicates that a compensation of the errors 
will, in a certain measure, he found in the values of the resultants 
of the actions. A solution of greater rigor would, however, he 
greatly desirable. 
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PART III. - DRAG AMD PROPULSIVE FORCE 12 

The component parallel to the velocity V of the aerodynamic 
action of the wing of infinite aspect ratio in nonsteady motion can 
"be readily computed on the basis of the solution of the problem of 
two-dimensional motion, as given in part I. Under the assumptions 
made, for any law of motion of the infinite wing, the theory per- 
mits computing the drag or propulsive force. For the hypothesis of 
a perfect fluid, the profile drag is not considered, nor are the 
variations of this drag due to the unsteady motion computable with 
the aid of this analysis. These actions are therefore added to 
those that are here computed. 

For the wing of finite aspect ratio, in the problem under 
condide'ration, the uncertainties mentioned in part II are also 
encountered. The analysis will therefore be limited to the results 
obtained for the two-dimensional motion. 1 ^ 

Drag and thrust in unsteady two-dimensional motion. - The 
symbol Rj is the instantaneous value of the force that arises in 
the direction of V on a segment of unit chord of the wing in a 
uniform flow of velocity V. With the notation of the preceding 
parts and in the same range of validity of the theory -there given, 

Rj., considered positive if it has the sign of a resistance and neg- 
ative if it corresponds to a propulsive force, can be computed with 
the aid of the egression given by Birnbaum (reference 3) 

R^ = PV^P ^ 7 dx - npa 2 L/4 (57) 

where 

2a = lim ^ 7 sin £ (58) 


12 

The numbers of the figures, the equations, and the paragraphs 
follow from the preceding part. 

1 " ) The treatment of Schme idler (reference 46) examines the aero- 
dynamic action corresponding to assigned vorticity of the wing. The 
method cannot, however, be generalized. 

14 The sign indicates the integral taken over the chord of the 

wing from the leading to the trailing edge. 
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The aerodynamic action on a segment dx of the chord is repre- 
sented by the force p Yy dx normal to the line of the axis and is 
therefore inclined to the 2 -axis normal to V by the angle dz/dx. 
The integral in the second member of equation (57 ) therefore repre- 
sents the action along x that is exerted on the points of the 
chord. The negative term in the expression corresponds in every 
case to a propulsive force and arises from the suction that is 
exerted at the leading edge by the surrounding fluid and produces 
a lowering in pressure that becomes infinite for the wing of infini- 
tesimal thickness* If 7 is expressed by means of the customary 
series of functions cot d/ 2 , sin d, . .., sin nd, only the first 
term can give rise to a suction at the leading edge because the 
other terms represent circulations that vanish. The quantity a 
defined by equation (58) gives the coefficient of the first term as 
found when considering that for d — > 0 , 


lim sin d cot jr = 2 


The steady motion that is considered corresponds to a chordwise 
distribution of vorticity represented by the same 7 as for the 
steady motion. Inasmuch as the suction at the leading edge depends 
on the instantaneous value of 7 , the value of 7 must be the 
same for nonsteady motion as for steady motion. In addition, when 
it is considered that for steady motion the resultant force in the 
wing direction must be zero. 


dx - flPa 2 L/4 (59) 


axis .of the wing on which, under 
the vorticity distribution 7 


where (dz/dx)o is the slope of the 
the considerations of steady motion, 
holds . 


If Wy denotes the velocity induced by these vortices, then 



w„ 
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(57) and (59), 


7 dz (60) 


which is a second form of equation (57 ) given by Jaeckel (refer- 
ences 28 and 29) and is entirely equivalent to the equation given 
in the development of the computations . 

Denoting by w \ the velocity induced by the free vortices and 
taking account of equation (1) yield 


From this relation and from equations 

El ' p Ji/ T s • w >') 


+ V! =y|| + || 

Equation (60), on the basis of this relation, becomes 


Hi 



W2 y dz - P 



dz 

St 


7 dz 


(61) 


This third form of equation (57) was used by Schmeidler (refer- 
ence 46). It lends itself to interesting interpretations: If the 

forces in the z -direction are distributed with density PV7 and 
the corresponding velocities of the points of application are Sz/St, 
the instantaneous power Nj absorbed by the motion of the points of 

the profile in the direction normal to V is ezpressed by the 
relation 



(62) 


It is noted that the values of z are assumed positive down- 
ward, whereas the pressures PV7 are positive upward. Hence Wj_, 

according to equation (62), is positive if work is done in overcoming 
the aerodynamic action. 
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When equation (62) Is considered. It Is observed that the 
second term in the second member of equation (61) represents the pro- 
pulsive force (or thrust) 15 that would arise if the phenomenon occurred 
without dissipation of energy; that is, without increase in the kinetic 
energy of the fluid surrounding the wing. The first term represents 
the drag R s that must he overcome by this phenomenon, or in other 
words, by the creation of the vortex wake. 

The computation of the instantaneous values Ri and Ni and, 
successively, of the mean values Ik and % of the same magnitude, 
in the case of harmonic motion, is immediately obtained on the basis 
of the expressions given in part I. Poggi (reference 47), on the 
basis of the investigations of Glauert (reference 20), computed the 
propulsion and power corresponding to a rotary oscillation and, as a 
limiting case, to the translatory oscillation of the wing. Ktissner 
(reference 21) determined these values for the wing in translational 
and rotary oscillation of the wing or flap. The same computations 
were made by Garrick (reference 48), who made use of equation (61) 
and of the energy interpretation of the tern 


R m = |pLV2 



J 

H-,/ 2 ) + iH 0 ( 2 ) 



(63a) 


- N m = IpLV 3 $> 


(63b) 


where 


15 This Is the quantity called "Vortrieb " by Schmeidler; the 
first term is called "Widerstand . " 
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$> = 



JJ X ^ 

l! (2) ♦ W 2 V 



lOt* 

and z = Ze is the ordinate of a point of the wing according to 
the oomplex notation with gt (C), |c{, and C the real part, the 
modulus, and the conjugate complex, respectively, of the complex 
quantity C, with H the Hankel function of reduced frequency 
U>=QL/2Y. 

The first addend that appears in the parentheses of the expression 
of % corresponds to the resistance Rg . If E denotes the ampli- 
tude of the sinusoid that represents the distribution of the free vor- 
ticity in the wake, it is found that, according to equation (63), the 
following relation holds: 


R s = pLe2/16 cq 


i " 

This relation can also he derived by considering that the work 
done by the resistance R s during the displacement 2jtV"/Q corres- 
ponding to an entire oscillation of the wing must be equal to the 
increment that the kinetic energy of the fluid has received in the 
same time, and therefore to the ♦kinetic energy of the fluid (considered 
stationary at infinity) in a strip included between two parallels to 
the z-axis at 2jtV/ Q distance from each other and located in the wake 
at a great distance downstream of the wing. 
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Appllcatlona to particular cases of oscillatory motion. In the 
case of a nondeformable wing (a vlng in motion of translation and 
rotation) , 


Z c Zm * Zq COS ^ 


in "which represents the complex amplitude of the oscillation of 
the middle point of the vlng and Zq the amplitude of the incident 
oscillation. Equation (63) yields, by simple computations. 


vhere 


M 2 ] 

(64a) 

M 2 ] 

(64b) 


J a Zo + iWZm + i iooZo 
2 


Hi + IHq 


H 1 (2) + iH 0 ^ 




+ iH 0 
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The quantities K, h, and k are functions of the reduced 
frequency. The variation of h and k as functions of u> are 
shown in figure 14. In figure 15, the real part is plotted on the 
abscissa of the quantities K/h and K/k and the imaginary part 
on the ordinate assuming the segment OA as unity and the positive 
sign of the imaginary axis in the downward direction. The points 
on the curve give values of the reduced frequency. 

In the particular case of translatory oscillation (Zq= 0), from 
equation (64), denoting by v »>^|Zm| the maximum velocity corres- 
ponding to this oscillation, 

-R m » 2. pLv 2 h (65a) 

N m = 2 pLv 2 kV (65b) 

The force that arises in this type of motion is a propulsive one. 
The efficiency of the wing considered as a means of propulsion is 
h/k, which is equal to 1 for u>=0 and decreases continuously toward 
0.5 as u> increases. 

If a wing of velocity vq normal to the wing velocity V is 

considered under conditions of steady motion, the lift P, the coeffi- 
cient of which is equal to nvq/V, under the usual assumptions gives 
a component in the dlfection Y represented by 

Pvj/V = jtpLv 2 1 


for a segment of unit chord. This component is directed forward 
whether vi is directed downward or upward. On varying harmo- 

nically, if it were valid to apply at each instant the expressions of 
the steady motion, vq would be the mean value of the propulsive force 
given by equation (65a), in which h = 1. It is therefore concluded 
that the exact analysis corrects this approximate consideration by 
reducing the propulsive force by a faotor depending on the reduced 
frequency. This factor approaches 1 when the reduoed frequency is 
decreased and approaches 1/4 when u> is decreased. The power absorbed 
is reduced according to a factor that varies from 1 to 0.5 with an 
Increase in ou . 



70 


NA.CA TM 1277 


It Is of interest to examine how, by combining a torsional motion 
with translational and oscillatory motion, the propulsive force can he 
increased. By simple computations, it is found, on the basis of equa- 
tion (64a), that for a given amplitude z p of oscillation of the rear 

neutral point of the wing, for every reduced oscillation frequency a 
certain value of the amplitude and of the phase of the torsional motion 
exist for which the propulsive force is a maximum. If K' and K" 
denote the real and imaginary parts of the quantity K/h, respectively, 
the component f of the rotation in phase with the translatory motion 
is given by 


f 

L K» -1 


(66a) 


The component in quadrature is expressed by 




Zp 2-K' 
L K'-l 


U) 


(66b) 


The propulsive force under these conditions is given by 


^m 


| pLv 2 h 


K ’ 2 + K " 2 
4(K'-1) 


(67) 


With respect to the propulsive force of the purely translatory 
osoillatlon expressed by equation (65a), the effect of the rotation 
introduces the factor dependent on the reduoed frequency 
(K ,2 +K" 2 )/4(K* -l) . This quantity assumes decreasing values with 
increases in the reduced 1 frequency, tending asymptotically to the 
value 1.125. For u) = 0.5, this value is equal to 1.445; for 
reduced frequencies not too small, the adding of the torsional motion 
does not greatly modify the value of the maximum obtainable propulsive 
force, the base value of which is always that of the purely flectional 
motion. For sufficiently small reduced frequencies, there are con- 
siderable increases. By neglecting the higher powers of the reduced 
frequency, the expression for % can be assumed 


*m«| ZpVV 
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On decreasing od, the torsional motion that must he combined with 
the flectlonal motion to obtain the maviminn propulsive force tends to 
assume a phase displacement of 90° ahead with respect to the trans- 
latory oscillation (that is, in the phase in which the translatory 
velocity is a maximum upward, the wing is nearest its ma yirmnn negative 
incidence angle) . 

The analysis of the variation of the propulsive force as a func- 
tion of the flectlonal motion for a given amplitude of the torsional 
motion was given in reference 49 by means of a graph that permits 
computing directly from the propulsive force (or drag) and the power 
absorbed (or emitted) in the oscillation. All the values of the 
efficiency from 1 to 0 can be obtained by suitably varying the ampli- 
tude and the phase of the flectlonal motion. The region of may i mum 
efficiencies is nearest the point J = 0, which corresponds to the 
motion without drag and without absorbed power, with zero vorticity 
in the wake. Hinder these conditions, for a small u>, the torsional 
motion is displaced by about a 90° lag with respect to the flec- 
tlonal motion; that is, in the phase in which the translatory velo- 
city^-® upward is a maximum, the wing is nearest the maximum positive 
incidence. This result, in relation to that of the analysis of the 
maximum propulsive force previously indicated, leads to the conclu- 
sion that the conditions of maximum efficiency are not compatible 
with those of maximum efficiency of the wing considered as a propul- 
sive means, which can be obtained with a certain loss in efficiency. 

On the basis of equation (63), it may also be determined whether, 
by combining with the flectlonal motion a deformation that alters the 
curvature, any advantage in the value of the thrust can be obtained. 

When 


then 


Z ** Zm + Z 2 cos 2d 


J = 2Zg + 






Hl (2) ^ 

h j 2) : w 8 v 



g PLV 2 


16 More precisely, the velocity of the rear neutral point. 



72 


NACA TM 1277 


In this case, results are obtained that are entirely analogous 
to those of the motion of translation and rotation. When 


2 Z 2 *= (f+iq) Znx 


there are obtained for f and q the same expressions of equa- 
tion (66) in which K' and K" are substituted for the real and 

imaginary parts of H^ 2 . With this modification, 
the factor of increase in the maximum thrust has the same expres- 
sion as for the preceding case. This factor, which can be expressed 
by l/4(k-h) , has the value for oj — >°° and the value 1.148 for 
w = 0.5. In this case also, for a not very small oo, no great 
increases are obtained in the thrust as compared with the purely 
translatory motion. As in the preceding case, for a given frequency 
of oscillation on increasing the wing velocity, then maximum thrust 
for a given translatory amplitude first increases rather slowly; 
only when sufficiently low values of the ratio Ql/ 2V are obtained 
does the maximum thrust tend to increase linearly with the velocity. 
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PART IV. - EXPERIMENTAL INVESTIGATIONS 

The experimental investigations that have thus far been conducted 
on the aerodynamic actions on a wing in unsteady motion are not as 
numerous as would he required by the complexity and Importance of the 
problem. It is from the measurement of the forces on the oscillating 
wing that conclusive data are expected that would permit a reliable 
computation of the critical velocities of the wings and tail surfaces. 
Various problems relative to the stresses of the wing structures 
during flight in agitated air also require experimental clarification. 
The experimental investigation should furnish the necessary control 
for the fundamental hypotheses of the theory of wings of infinite 
aspect ratio and for the finite wing, the actual theory that makes 
vise of approximations that have not yet been completely checked should 
be integrated. The research presents, in addition to the difficulties 
common to all problems for which forces variable in time are to be 
measured, serious obstacles for the requirement of absolute regularity 
of the stream in which the experiment is conducted. Small fluctuations 
in the velocity and in the direction of the wing, which do not have 
any great effect in the measurements of a steady flow, can render the 
measurements of the forces on the oscillating wing entirely unreliable. 

In this part, the results obtained up to the present by various 
experimentors will be discussed, and the results compared with theory. 

English tests. - The first series of tests was. conducted by 
Duncan at the National Physical Laboratory and published in 1928 
(reference 50) . The object of the tests was to check the mechanical 
theory of the wing oscillations. From these measurements Duncan 
obtained, for a particular wing model, the values of the aerodynamic 
coefficients to be introduced in the expression of the velocity in 
order to compare the calculated value with the experimental velocity 
obtained. The greater part of the tests was conducted on a model 
that was deformed during the oscillation, according to an incompletely 
defined law. The tests therefore do not lend themselves to a check 
of the aerodynamic theory, a check with which the experimentor was 
not concerned, as he did not then have the results of the theory. 

A series of tests were, however, conducted by Duncan on a model that, 
during the oscillation, rotated rigidly about an axis parallel to the 
span. The wing was rectangular, with RAF 15 profile, 152-millimeter 
chord, and 686-millimeter span. The axis of rotation was at l/lO of 
the chord from the leading edge. The damping of the oscillations was 
measured in the presence of a wind and in still air for various angles 
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of attack and frequencies of oscillation. Inasmuch as the oscilla- 
tions did not have very rapid damping, the results can be compared 
with those of the theory on harmonic motion. According to the theory, 
the moment of the aerodynamic force due to the rotational oscillation 
possesses a component in phase with the motion and a component Mq 

in quadrature and therefore in phase with the angular velocity q; 
the moment Mq, which has its sign opposite to q, therefore consti- 
tutes a damping action and may be put in the form 


Mq » Jtpbq L 2 VS 

where S is the wing area. For a segment of an infinite wing, the 
coefficient b depends on the reduced frequency and on the position 
of the axis of oscillation. On increasing the reduced frequency, b 
tends to the value 17 


( 68 ) 


where t la the distance of the axis of rotation from the focus con- 
sidered positive if the axis is in the rear. For the finite wing, b 
depends also on the plan form. Its values for not-too-small reduced 
frequencies are not, however, considerably removed from that given by 
equation (68). For the position of the axis of rotation of the tests 
of Duncan, this gives b = 0.211. The values of b obtained on the 
basis of the damping moments measured by Duncan are given in figure 16. 
For the computation of b, the value of the aerodynamic damping 
moment is considered to be the difference between the measured values 
in the presence of wind and in still air. The values are all below 
that given by equation (68). For equal velocities, these values Indi- 
cate a decrease with a decreasing u>, as would also be given by theory 
for positions of the axis of rotation ahead of the focus. For equal 
frequency, on decreasing the velocity (hence on increasing co), the 
experimental values in general, indicate a decrease that can be 
ascribed to the effects of the Reynolds number. There is also a 
decrease on decreasing the angle of attack (at least in the region 
investigated). For angles of attack from -4° to -5°, Duncan found 


17 Which can be derived on the basis of the expressions given in 
references 2 and 21. 
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a vanishing of the aerodynamic damping. This phenomenon, called by 
Studer "oscillations of separation," cannot be studied by the theory 
of the preceding parts. 

The measurements conducted by Duncan of the damping of the 
oscillations of the flap also lend themselves to a comparison with 
theory. In these tests, the damping due to the friction of the 
suspension is rather large. By assuming for this case also that the 
aerodynamic damping can be obtained from the difference between that 
measured with wind and that measured in still air, it is found that 
the experimental value of this damping is equal to about one •half the 
theoretical. This disagreement should not be surprising, because the 
derivatives relative to the flap are always markedly less than the 
theoretical values. In the tests conducted by Duncan, the value of 
the derivative of the hinge with respect to the angle of the elevator 
under conditions of steady motion was equal to 0.6 of the theoretical 
value . 


Tests conducted at the Laboratorlo dl Aeronautloa dl Torino. - 
The aerodynamic actions on the oscillating wing were measured by Cicala 
in the free -jet wind tunnel of 600-millimeter diameter at the 
Laboratoria di Aeronaut ica di Torino. The chord of the models on 
which the tests were conducted was about 13 centimeters and the span 
about 50 centimeters. Because of the relatively small dimensions of 
the jet, which was free in the region in which the model was located, 
the wing operated with a rather low effective aspect ratio: the 
value of SCp/da under steady conditions (referred to pV^) was equal 

to about 0.5« because the wing projected from a plane that was placed 
tangent to the Jet in order to mask the suspension and measuring 
apparatus. This rather low value of the aspect ratio is one dis- 
advantage of these tests, which are described in references 14, 51, 
and 52. 

In a first series of tests (reference 14), the aerodynamic 
damping of the flectional oscillations (rotation of the model about 
an end chord) and the damping of the oscillations about the axis con- 
taining the foci of the various sections (also a rigid rotation) were 
measured. The measurements relative to the flectional motion were of 
little importance because of imperfections in the construction of the 
model and of the measuring apparatus. These tests were later repeated 
(reference 52). The measurements relative to the torsional oscilla- 
tions gave for b a value of about 0.11 (against 0.125 given by equa- 
tion (68) for |=0), which was almost constant when the reduced 
frequency was varied as required by theory. This value was confirmed, 
at least for the range of not very large angles of attack, by succeeding 
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tests on a model different from the one described in reference 14, 
this model also being of symmetric profile but of greater rigidity. 

The principle of the measuring apparatus for these tests and for 
those of series III can be briefly described as follows : The oscil- 

lation of the model was controlled through an intermediary element 
that possessed two simultaneous motions, a rotation depending on the 
displacement imposed on the model and a rotation about an axis per- 
pendicular to the first rotation and depending on the magnitude of 
the force transmitted. This element carried a mirror that reflected 
on sensitive paper a luminous point that, by describing the motion 
of the intermediary, gave the force -displacement diagram (and there- 
fore the moment-rotation diagram). The test conducted for equal fre- 
quency in still air and in the presence of wind permitted isolating 
the aerodynamic action. In figure 17 are given some of the oscillo- 
grams thus obtained that give the simultaneous values of the angular 
position of the wing and of the moment transmitted. The field of 
the coordinates can be retained as Cartesian, so that the diagrams 
are approximately ellipses. The enclosed area measures the work 
absorbed in the oscillation. In the figure are given the scale of 
motions and also the lines a = constant corresponding to the extreme 
positions for one of the Oscillograms obtained for a wing velocity of 
9.4 meters per second and for a number of oscillations eq ua l to 570. 

For each measurement, the oscillogram was obtained by permitting the 
luminous point to run through two or three cycles. The paper was 
successively advanced, thus intercepting the light point. In order 
to obtain a reference point of the angles of attack for each oscillo- 
gram, a point in a fixed position was marked (points F in the figure). 

In the third series of tests (reference 52), the damping measure- 
ments of the translatory flectional oscillation (in the sense pre- 
viously defined) were repeated. The component of the lift in phase 
with the displacement of the wing can be expressed by the derivative 


1 d g p _ 

n <3(v/Y) 2 


where v is the translational velocity normal to V and v/v is 
the corresponding variation of the angle of attack. According to the 
theory of the infinite wing, the factor ag varies from -1 to -0.5 
on increasing the reduced frequency (always negative because the vari- 
ation of the lift has a sign opposite to the vertical velocity). 
According to the approximate theory of the finite oscillating wing, 
the range of variation of this factor is reduced, starting from 
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^ = 0, from the corresponding value of the steady motion computed 
with the aid of the vortex -filament theory of Prandtl. The tests 
gave for ag a value nearly constant and equal to 0.43 for not 
too small a velocity within the range of reduced frequency in which 
the tests were conducted (0.2 < oj <0.6). At small velocities, a 
decrease of ag was found that was ascribed to the effect of the 

low Reynolds number. 

The lift component in phase with the fleet ional motion and the 
focal moment in phase with the torsional motion were rather small, 
as required by theory, the aerodynamic inertia effect being included 
with the measurements of the mass of the model on the basis of the 
oscillation data in still air. 

The variations of lift that arise from an oscillation of the 
wing about the focus were also measured (reference 51) . By expressing 
the lift component in phase with the rotation by means of the derivative 



values of about 0.52 were obtained for this coefficient for all the 
reduced frequencies at which the tests were conducted (0.1 < w < 0.7). 

The component Pq of the lift, proportional to the angular velo- 
city q and therefore in quadrature with the motion, can be expressed 
by means of the relation 

Pq e> «pa4qSLV 

The coefficient a 4 , which, according to the theory of the wing 

of infinite aspect ratio, has values increasing with u> and approaching 
to 0.5, was found from the tests to be almost always equal to 0.38. 

The principle of the apparatus for the measurement of the lift 
due to the torsional moment was the following: The wing was put in 

forced torsional oscillation by guiding, according to the harmonic 
law, a point of an end section while a point of the same section was 
attached by means of a steel wire. Under these conditions, in addi- 
tion to the torsional motion, a flexional motion of rotation arose 
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about the chord of this section. In another section, a force having a 
component in phase and one in quadrature with the excited motion, was 
introduced, the amplitudes of which could be varied during the test. 

In the presence of wind, the amplitudes of the two components were 
controlled so as to eliminate the flectional motion and to balance the 
aerodynamic action and the inertia force. From the tests conducted at 
equal frequency in still air and with wind, the aerodynamic forces were 
obtained. By the same principle, the focal moment due to the flectional 
motion was measured. In agreement with theory, the value of the 
flectional motion was so small that it could not be measured . 

For all the aerodynamic derivatives expressed in the preceding 
form, almost constant values were thus obtained in these tests by 
varying the reduced frequency. The theory of the wing of finite aspect 
ratio developed in reference 37 justifies this result for the range 
of not-too-small reduced frequencies and low aspect ratios at which 
the tests were conducted. In figure 16, the dotted line gives the 
value of the coefficient b that would be obtained by these tests, 
a value that does not diverge much from the tests of Duncan. 

American testa. - Tests have recently been conducted in the 
United States for the measurement of the aerodynamic forces on the 
oscillating wing to check the theory, of the infinite wing. An 
interesting series of tests was conducted by Reid and Vincenti at 
the Guggenheim Laboratory (reference 53). The model used had a chord 
of 38 centimeters and therefore permitted the attainment of suffi- 
ciently high Reynolds numbers. The span was not large (about 91 cm). 
Nevertheless, a large aspect ratio was obtained because the model was 
placed between two walls normal to the plane of the wing. The wing 
of NACA 0015 profile was put in oscillation about an axis at a dis- 
tance of 4/10 chord from the leading edge. At the opposite edge to 
that at which the motion was excited, the aerodynamic action was 
measured. The wing support, consisting of a ball bearing, was sus- 
tained by a rigid spring the inflections of which were recorded by 
means of mechanical and optical amplification on a strip of sensi- 
tive paper with uniform forward motion. On the same strip were marked 
the instants at which the wing occupied the extreme and middle posi- 
tion. With the aid of a harmonic analysis of a graph of the forces, 
which were necessarily irregular, the amplitude and the phase of the 
fundamental harmonic with respect to the motion of the wing were 
derived and thus the ratio r of the amplitude of the lift under 
conditions of oscillatory motion and the ratio corresponding to the 
steady motion for equal rotation and phase angle (leading) 5 
between the lift and the rotational motion were measured. The 
results are plotted as a function of the reduced frequency in 
figures 18 and 19 and compared with the theory of wings of infinite 
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aspect rations (continuous curve) and with those that were obtained 
in the Torino tests. The Torino tests must necessarily present a con- 
siderable divergence from the American tests because of the difference 
in aspect ratio. The phase displacements predicted by theory are 
somewhat greater than the experimental values. A considerable devi- 
ation is presented by the theoretical and experimental curves of the 
ratio r. 

A similar series of tests was conducted by Silverstein and 
Joyner (reference 54). The model had a chord of 13 centimeters, 
considerably smaller than that used in the tests by Held. In this 
case also, a large aspect ratio was attained by using end walls. 

The oscillation axis passed through the forward quarter chord and 
the aerodynamic lift force was measured by means of an apparatus 
based on the same principle as previously described. In these tests, 
only the phase displacement between the lift and the rotation was 
measured and the values shown in figure 20 were obtained; the con- 
tinuous lines give the values of the theory of the infinite wing and 
the dotted line give the values obtained from the Torino tests. The 
scatter of the test points is large for the high values of w ; that 
is, for the tests conducted at low velocity. 

Because of the small number of the results that are available, 
no deductions of a conclusive character can be given. The different 
conditions under which the tests were conducted also does not provide 
a good basis for comparing the different results. 

The tests in which the conditions for a check of the theory of 
the oscillating wing of infinite aspect ratio were best realized are 
those of the Guggenheim laboratory. The comparison is not, however, 
completely satisfying. The probable cause of the divergence encountered 
seems to lie in the agglomeration and dissipation of the wake vortices, 
the mutual positions and intensities of which the theory assumes to 
be maintained indefinitely. The problem should be investigated more 


18The theoretical curves of the graphs of reference 53 do not 
coincide with those given in figures 18 and 19 because the aerodynamic 

2 

inertia effect represented by the terms in u) in the expressions of 
the derivatives is not considered. In fact, this action, which remains 
unchanged with and without wind, is already compensated in the pre- 
liminary operation of putting the center of gravity on the axis of 
rotation, a compensation that, it seems, was effected under dynamic 
conditions. The correction is small and makes the test points 
approach the theory more closely. 
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thoroughly, especially in an experimental manner. The verification 
of the theory of the infinite wing is less urgent, however, th^ n the 
investigation of the finite wing, particularly for the phenomenon of 
wing vibration, in which, because the motion is more pronounced toward 
the tip of the wing, the conditions are considerably removed from two- 
dimensional motion. This limiting case is also difficult to obtain 
experimentally because of the considerable importance assumed by the 
wake over a large distance behind the wing. The Torino tests make 
use, however, of low aspect ratios for which the approximations of 
the theory of reference 31 are less justified for giving an account 
of the results of such tests . There would therefore be- required : 

First, a perfecting of the theory of the oscillating wing of finite 
aspect ratio; and second, the extension of tests to wings of greater 
aspect ratio. The range of angle of attack within which the coeffi- 
cients can be held constant must be defined and the field of coeffi- 
cients relative to the oscillating wing with flap must be investigated. 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics. 
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TABLE I 


a 

R 

a 

R 

a 


0 

0.6000 

1.7 

0.3490 

7.6 

0. 1688 

0.1 

0.4878 

1.8 

0.3427 

8.0 

0. 1609 

0.2 

0.4762 

1.9 

0.3366 

8.6 

0. 1436 

0.3 

0.4661 

2.0 

0.3307 

9.0 

0.1368 

0.4 

0.4646 

2.1 

0.3250 

9.6 

0.1307 

0,6 

0.4443 

2.2 

0.3196 

10 

0. 1250 

0.0 

0; 4346 

2.3 

0.3141 

11 

0. 1147 

0.7 

0.4253 

2.4 

0.3088 

12 

0.1058 

0.8 

0.4163 

2.6 

0.3038 

16 

0.0862 

0.9 

0.4077 

3.0 

0.2806 

20 

0.0634 

1.0 

0.3994 

3.6 

0,2600 

25 

0.0499 

1.1 

0,3914 

4.0 

0.2420 

30 

0.0408 

1.2 

0.3837 

4.6 

0.2261 

40 

0.0298 

1.3 

0.3763 

6.0 

0.2118 

50 

0.0232 

1.4 

0.3691 

6.6 

0.1990 

100 

0.0109 

1.5 

0.3622 

6.0 

0.1875 

600 

0.0020 

1.6 

0.3666 

6.6 

0. 1770 

1000 

0.0010 



7.0 

0.1676 

00 

0.0000 
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TABLE II 


R, 

0.0000 
0.2824 
0. 3961 
0,4812 
0.6613 
0. 6116 
0.6649 
0.7128 
0.7563 
0.7964 
0.8334 
0. 8678 
0. 9001 
0.9303 
0.9588 
0.9858 
1.0113 



s, B, 


8,5 

1.6642 

9.0 

' 1.6817 

9.5 

1.6976 

10 

1.7123 

10,5 

1 

12 

1.7603 

16 

1.8235 

20 

1.8624 

25 

1.8934 

30 

1.9135 

40 

1.9380 

50 

100 


500 

1.9958 

1000 

1.9980 

00 

2.0000 
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TABLE I 


o > 

i — x 

X‘ 

0 

1,0000 

0 

0.002 

0.9967 

0.0126 

0.01 

0.9826 

0.0456 

0.02 

0.9637 

0.0752 

0.04 

0,9267 

0.1160 

0.06 

0, 8920 

0. 1426 

0.08 

0,8604 

0.1604 

0.10 

0. 8319 

0.1723 

0.12 

0. 8063 

0. 1801 

0. 14 

0. 7834 

0.1849 

0.16 

0. 7628 

0.1876 

0,18 

0. 7443 

0,1887 

0.20 

0. 7276 

0.1886 

0.22 

0.7125 

0,1877 

0.24 

0.6989 

0. 1862 

0.26 

0. 6865 

0.1842 

0.28 

0. 6752 

0.1819 

0. 30 

0.6650 

0,1793 

0. 32 

0. 6556 

0.1766 

0.34 

0, 6469 

0.1738 

0,36 

0.6309 

0.1709 

0, 38 

0.6317 

0.1679 

0,40 

0,6250 

0.1650 

0.42 

0.6187 

0.162! 

0.44 

0,6130 

0. 1592 

0.46 

0.6076 

0, 1563 

0.48 

0,6026 

0.1535 

0.50 

0.5979 

0.1507 

0,52 

0.5936 

0. 1480 

0.54 

0.5895 

0. 1454 

0.66 

0.6857 

0. 1428 

0,68 

0.5822 

0. 1402 

0.60 

0.5788 

0.1378 




U) 

1 — X 

X' 

0,62 

0,5757 

0. 1354 

0.64 

0. 5727 

0. 1330 

0.66 

0, 5699 

0.1308 

0,68 

0,5673 

0. 1286 

0,70 

0,5648 

0.1264 

0.72 

i 

0.5624 

0.1243 

0,74 

0, 5602 

0,1223 

0. 76 

0, 5681 

0.1203 

0. 78 

0. 5561 

0,1184 

0.80 

0,5541 

0.1165 

0, 82 

0,6523 

0,1147 

0.86 

0, 5490 

0.1112 

0,90 

0. 5459 

0, 1078 

0.94 

0.5432 

0.1047 

0, 98 

0.5406 

0, 1017 

1.00 

0, 5394 

0,1003 

1.1 

0, 5342 

0,0936 

1.2 

0. 5300 

0.0877 

1.3 

0,5265 

0,0825 

1.4 

0.6235 

0.0778 

1.5 

0, 52 JO 

0,0735 

1.6 

0.5 J 89 

0,0697 

1.7 

0.5171 

0,0663 

1.8 

0,5165 

0.0632 

1.9 

0,5142 

0.0603 

2,0 

0,5130 

0.0577 

2.6 

0, 5087 

0.0473 

3.0 

0,5063 

0.0400 

3.5 

0,6047 

0.0346 

4.0 

0,5037 

0.0306 

4,5 

0, 5029 

0.0273 

5,0 

0.5024 

0,0246 

10,0 

0. 5006 

0,0124 

>10 

0. 5000 + 

+ 1/16S* 

l /8« 
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